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L  INTRODUCTION. 

TUf  report  preeente  some  nanrigation  and  almanac  algorithma  with  implementation  pro* 
grami  that  are  written  in  the  Radio  Shack  TS8-ao  MODEL  4  BASIC  language.  This  language 
can  be  customised  for  the  Nacry  adopted  Sharp  PO-ISOOA  or  other  BASIC  language  com¬ 
puters. 

Because  the  PC-isooA  has  limited  memory,  many  of  the  classical  narigation  formulas 
horn  spherical  trigonometiy  hare  been  rewritten  in  a  fonn  which  minimises  or  eliminates 
the  need  to  determine  special  cases  especially  those  of  quadrant.  Det^  are  contained 
in  Section  IL 

The  navigation  algorithms  are  detailed  in  the  BASIC  program  NAVAL60R,  presented  in 
Section  DL  Programs  for  the  determination  of  the  position  of  the  Sun,  Moon,  Venus,  Mars, 
Jupiter  and  Saturn  are  contained  in  the  BASIC  program  NAVEPHM  which  is  presented  in 
Section  IV.  These  programs  reproduce  the  tables  in  The  Maieiical  Almanac  and  The  Air 
AlmoMac  to  within  0.2^. 

n.  THE  GENERAL  SPHERICAL  TRIANGLE 
AND  COORDINATE  TRANSFORMATIONS. 

A.  Backgroimd.  Many  of  the  algorithms  and  procedures  used  for  navigational  compu¬ 
tations  were  designed  for  use  with  tables  of  tines,  cosines  and  tangents.  To  save  space  and 
diminate  repethic!!,  these  tables  usually  contained  only  values  for  the  first  quadrant.  To 
generate  values  for  other  quadrants,  many  cumbersome  rules  were  adopted.  These  rules — 
contained  in  classical  references  such  at  Bowditch  [Ref.  1] — are  promulgated  in  modem 
computer  programs,  causing  needlessly  additional  programming  effort.  This  additional  ef¬ 
fort  can  largdy  be  eliminated  by  using  the  concept  of  the  general  spherical  trian^^e  whidi 
is  described  in  the  next  section. 

B.  Theory.  The  formulas  for  the  general  spherical  triangle  expounded  on  by  Wm.  Chau- 
venet  [Ref.  3]  in  his  book,  A  Tnetin  on  Flans  and  Spherical  Trigonomeirif,  which  was  first 
published  m  1850.  One  usually  thinks  of  a  spherical  triangle  as  looking  something  like  the 
object  depicted  in  Figure  1,  where  the  arc  loigth  of  each  leg  is  leas  than  180*.  Chanvenet’s 
general  triangle  does  not  have  the  180*  limitation  and  so  objects  nimilar  to  the  one  shown 
in  Figure  2  are  also  considered  to  be  spherical  triangles.  In  Figure  2.  the  croespoint  P 
is  180*  from  the  vertex  of  angle  A  In  labeling  spherical  triangles  it  is  important  that  an 
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angle  be  labeled  whh  a  capital  letter  and  that  the  opposite  aide  be  labeled  with  the 
bnt  lower  case  letter. 

Depending  upon  the  quantitiea  to  be  determined,  either  the  set  of  equations 

cos  a  s  cos  i  cos  c  +  sin  6  sin  c  cos  ii, 
sinacos^  cosisine  -8in6cosccoSi4, 
sin  a  sin  B  s  sin  6  sin  ii, 

or  the  polar  form  of  Equs.  (1) 

cosii  =  -  cos  B  cos  C  +  sin  B  sin  (7  cos  a, 

8inilcos6  =:  cosBsinC  +  sinBcosCcosa, 


nuqr  be  need. 

ChamreBet  stated  that,  *It  will  be  fomul  that  aUtk»  tu  eases  of  the  general  triangle 
edmd  two  eolniione,  bat  that  tkef  all  become  ietermmaie,  when,  in  adUHon  to  the  other 
data,  the  sign  of  the  tine  or  eoeine  of  one  of  the  parte  it  given.* 

fbnniilaa  to  spherical  co(»dmates  shonld  be  developed  so  that  quantities  snch  as 
latitude,  declmatioii  and  altitude,  which  have  vahiea  in  the  range  -00*  to  +90*,  are  de* 
tennined  bjr  arcame  or  arctangent  fommlas;  so  that  quantities  such  as  longitude,  aainmth, 
Greenwich  hour  an^  and  local  hour  an^  whidi  have  values  in  the  range  -180*  to  +180* 
or  0*  to  360*,  are  determined  by  both  a  sine  and  a  cosine  formula  used  in  conjunction  with 
a  quadrant  detennining  arctangent  function  (the  qatn  function);  and  so  that  a  quantity 
snch  as  distance,  which  has  a  spherical  arc  value  in  the  range  of  0*  to  180*  (assuming  the 
user  wishes  the  shortest  distance),  is  determined  by  an  arccosine  formula. 

Hhistmtkm  1.  Consider  the  *inverBe*  problem  of  spherical  geometry,  which  is:  Given 
the  latitude,  and  longitude.  A,  of  two  points  /^(^i.Ai)  and  determine  the 

distance  d,  the  forward  asimnth  au  and  the  badnrard  asnnnth  ot^i.  Tlie  solutions  will  be 
derived  using  the  convention  that  eastern  kmptndes  are  positive  and  western  longitudes 
are  negative.  Also,  northern  latitudes  are  positive  and  southern  latitudes  are  negative. 
The  geometry  is  iUnstrated  in  figure  3. 
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Sohre  for  d  and  013.  Label  the  anglea  and  legs  so  that  A  —  -  Xt^  B  —  013, 

C  =  960*  -  031,  a~d,6  =  90*-^  and  e  =  90*  -  ^1.  Sabetitnting  theae  into  Eqtu.  (1) 
and  limplifying  the  differencee  of  an^^,  we  obtain 

cosdssin^an^i  +coe^coe^i  coe(A3  -  At), 

8indcoeai3  sain^coa^i  -cos^sin^i coe(A3  -  At), 
sin  dain  013  =  coe  ^  ain(A3  -  At). 

The  first  equation  can  be  used  to  sohre  for  d.  Using  the  piinciple  angle  of  the  arccosine,  d 
lies  in  the  range  0*  to  180*.  If  the  non-prindple  an^e  is  used,  d  lies  in  the  range  180*  to 
960*,  we  can  either  travel  the  shortest  great  aide  distance  from  Pi  to  P2  or  we  can  travel 
the  great  circle  ronte  which  goes  around  the  backside  of  the  earth — the  fonner  solution  is 
assumed  to  be  preferred.  Hence,  d  is  restricted  to  lie  between  0*  and  180*.  Multiply  d  by 
60  to  get  the  distance  in  nautical  miles. 

Chanvenet  states  that  this  system  of  equations  has  two  solutions,  but  has  only  one 
solution  if  the  sign  of  the  sine  or  cosine  of  one  of  the  parts  is  known.  Since  d  lies  between 
0*  and  180*,  the  sign  of  sin  d  is  known,  thus  the  system  has  only  one  solution.  Farther, 
and  most  important,  the  sign  of  sind  is  positive  for  all  d  between  0*  and  180*  and  so  there 
are  no  q>ecial  cases. 

The  first  equation  is  used  to  compute  d: 

dsarccosjsin^sin^t +cos^cos^icos(A3~Ai)].  (9) 

The  second  and  third  equations  are  used  to  compute  013.  A  natural  tendency  would  be  to 
solve  either  the  second  equation  or  the  third  equation,  but  doing  so  loses  the  information 
concerning  the  quadrant  of  013.  Use  the  second  equation  to  determine  cos  013  and  the  third 
equation  to  determine  sin  013.  With  both  the  one  and  the  cosine  known,  the  quadrant  is 
uniquely  determined  using  the  quadrant  determining  arctangent  function.  The  solution  is 


013  =  qatn(sm  013,  cos  013) 

as  long  as  d  is  not  0*  or  360*.  Even  though  the  chance  of  division  by  sero  with  the 
occurrence  of  a  value  for  d  of  exactly  0*  or  exactly  360*  b  very  remote,  it  can  be  elimmated 
entirdy.  Since  the  sign  of  sin  d  is  positive  or  sero  in  both  the  second  and  third  equations, 
it  effectively  cancels  out  usmg  the  qatn  function.  The  asimuth  013  is  found  solving 


013  =  qatn(cos^8in(A3  -  Ai),8in(^coe^i  -cos^sin^i  cos(A3  -  Ai)].  (4) 


lb  find  the  back  aBnmth,  labd  the  angles  and  legs  ao  that  A  »  A3-A1,  B 
C  =  dta,  a=sdyb  s  90*-^i  and  e  s  90*-^.  Sabstitnte  these  into  Eqns.  (1)  and  simplify 
to  obtain 

asi  =  qatn[-cos^i  8in(Ai  -  A3),8in^i cos^  - coe^  sin^ cos(Ai  -  Aa)].  (5) 

mnstratioa  3.  Consider  the  ‘direct*  problem  of  spherical  geometry,  which  is:  Given 
the  latitode,  and  longitiide,  A,  of  a  point  FiC^iiAi),  the  distance  d  and  the  forward 
asimnth  aia  to  a  point  Aa),  determine  Aa  and  the  backward  asimnth  oai. 

Solve  for  ^  and  Aa.  Eelabd  the  an^es  and  lep  so  that  A  s  aia,  B  ^  Xi  -  Aa, 
C  =:  3fi0*  -  ttiai  asQO*-^,  and  c  —  00*  -  ^1.  Snbstitating  these  into  Eqns.  (1) 
and  simplifying  the  differences  of  angles,  we  obtain 

sm^  s  cosdsin^i  +sindcos^icosaia, 
cos^oos(Aa  -  Ai)  =  cosdcos^i  -sindsin^)  coeaia, 
cos^sm(Aa  -  Ai)  ^sindsinaia. 

The  first  equation  can  be  used  to  solve  for  Using  the  principle  angle  of  the  arcsine,  ^ 
lies  in  the  ranp  -00*  to  00*,  which  is  correct  for  a  latitude.  Smce  ^  lies  between  -90* 
and  +00*,  the  sign  of  cos  ^  is  known,  thus  the  system  has  only  one  eolation.  Farther,  and 
again  most  important,  the  rign  of  cos^  is  positive  for  all  ^  between  -90*  and  +90*  and 
so  there  are  no  special  cases. 

The  first  equation  is  used  to  compute 

^  =  arc8in(coedan^i +8indcoe^i  cosoia).  (6) 

The  second  and  third  equations  determine  coe(Aa  -  Ai)  and  8in(Aa  -  At).  As  before,  with 
both  the  sine  and  the  cosine  known,  the  quadrant  is  uniquely  determined  using  the  qatn 
function.  Since  cos  ^  is  always  positive,  it  can  be  eliminated  in  both  equations.  Then 

Aa^At  +qatn(8indsinaia,co6(fcoe^i  -  sindsin^i  cosoia).  (7) 

The  back  asimnth  can  be  found  by  relabeling  the  angles  and  legs  so  that  A  =  oai , 
B  =  360*  -  Oat,  C  =  At  -  Aa,  u  =  90*  —  ^,  h  =  90*  -  and  e  =  d.  Substituting  these 
into  Equs.  (1)  and  simplifying  the  differences  of  angles,  we  obtain 


Oat  =qatn(-cos^t  sinataisin^isind-cos^i  cosdcosoia). 


(8) 


Tbis  ncond  fommla  for  the  back  aihmith  ia  naed  in  the  spirit  that  the  enknoirn  parameters 
shonld  be  determined  aa  a  function  of  the  known  panunetera  rather  than  determining 
the  second  unknown  parameter  as  a  function  of  the  pierionaty  evaluated  first  unknown 
parameter.  That  is,  any  one  of  the  unknowns  can  be  found  as  a  function  of  the  knowns 
without  having  to  first  find  one  of  the  other  unknowns. 

mnstratkm  S.  Although  derived  tern  different  principles,  the  equations  ai  coor¬ 
dinate  transformation  obey  the  rule  of  Ghauvenet.  fbr  example,  the  equations  whidi 
transform  local  hour  angle  and  dechnation  to  altitude  and  asinmth  are 

cosasinA  »  -  cos ^rinh, 

cos  a  cos  A  sin  j  cos  ^  -  cos  £  sin  ^  cos  A,  (9) 

sina  s  sin  5  sin  ^ + cos  ^  cos  ^  cos  A, 

where  a  s  altitude,  A  =  asinmth,  S  »  declination.  As  local  hour  angle  and  ^  =  latitude 
[Ref.  4,  pg.  26}. 

The  third  equation  is  used  to  solve  for  altitude 

as  arc8m(sinj8in^+cos^cos^cosA).  (10) 

The  altitude  lies  between  -90*  and  +90*,  and  so  is  correctly  determined  by  the  arcsin 
function.  Since  cos  a  is  always  positive  it  can  be  eliminated  between  the  first  two  equations 
when  timng  the  qatn  function.  The  local  hour  an^e,  which  lies  in  the  range  0*  to  360*,  is 
determined  by 

A  =  qatn(- cosasinA,  sin  5  cos  ^ -cos  5  sin  ^  cos  A).  (11) 

C.  SfanpMeation  of  Standard  Works.  Using  the  general  triangle  and  the  qatn  func¬ 
tion,  the  simplifications  listed  below  can  be  made  in  Volume  II  of  Bowditch  [Ref.  Ij.  The 
various  sections  are  preceded  by  the  section  symbol  §  and  equation  numbers  from  Bowditch 
within  those  sections  are  prefixed  with  the  letter  B.  Unprefixed  equation  numbers  are  those 
contained  in  this  document. 

§706.  Sohring  for  altitiide. — Using  the  conventions  that  northern  latitudes  and 
declinations  are  positive  and  that  southern  latitudes  and  decimations  are  negative,  the 
altitude  of  a  star  can  be  computed  directly  from  Equ.  B(2a)  or  Equ.  B(2b).  Equ.  B(2a)  is 
equivalent  to  Equ.  (10).  All  of  the  special  cases  are  thus  eliminated. 

§707.  Sohing  for  ssiiimth.— Instead  of  using  Equ.  B(4b)  or  Equ.  B(5b),  use 
Equ.  (11).  See  §708  below. 


$TOS.  TIdm  asimatk* — Bowditch  oiet  meridiaa  an^es  which  can  be  labeled  either 
eaft  or  west.  Since  aititnde/aiinmth  is  a  left-handed  coordinate  system,  the  contention  is 
that  west  meridiaa  an^es  are  poothre  and  east  meridian  angles  are  negatite. 

BumpU  i.— The  latitude  of  the  obsoter  is  30*25/0N;  the  declination  of  the  celestial 
body  is  22*06/2N;  the  local  hour  angle  is  30*54/7W.  Using  Equ.  (11), 

A  -  qatn(-cos2r06'28in38*54/7,ain22*06/2cos30*25.'0 

-  coo  22*06/2  sin  30*25/0 coo  39*54/7) 

=  qatn(-{0.92661)(0.64161),  (0.37628)(0.86237) 

-  (0.92651)(0.50628)(0.76703)] 

=  qatn(-0.59446, 0.32449  -  0.35980) 
s  qatn(-0.59445,  -0.03531) 

-93.*39913 
»286.*60087 
A  =  266*36/1. 

BxampU  6.— The  latitude  of  the  observer  is  30*2S/0S;  the  declination  of  the  celestial 
body  is  22*06/2N;  the  local  hour  an^  is  39*54.'7K  Using  Sqn.  (11), 

A  =  qatn(-coo22*06/28in(-39*54/7),ain  22*06/2  coo(-30*25/0) 

-  cos  22*06/2  8in(-30*25/0)  coe(-39*54/7)l 
--  qatn(-(0.92651)(-0.64161),  (0.37628)  (0.86237) 

-  (0.92651)(-0.50628)(0.76708)1 
=  qatn|0.50445, 0.32449  -  (-0.35980)] 

=  qatn(0.59445, 0.68429) 

=  40.*98141 
A  =  40*58/9. 

D.  Spheroid  Earth  and  Great  Cbde  Fomnilaa. 

1.  Spheroid  Earth  Fomnilas.  The  spheroid  earth  algorithms  are  those  of  Thomas 
[Ref.  6|.  Some  changes  by  Shudde  [Ref.  6}  make  quadrant  determination  automatic  by 
use  of  the  qatn  function.  East  longitudes  and  south  latitudes  are  negative.  The  earth’s 
equatorial  radius  is  denoted  by  a«  and  the  earth’s  flattening  factor  is  denoted  by  /. 


a.  ImroM  Sohttkm.  Input  variables  are  the  latitude  and  lon^tude  Ai  of  Pi, 
the  latitude  ^  and  longitude  As  of  Pt,  ^  and  As.  Output  variables  are  distance,  (5/a«), 
forward  aiiznnth  from  Pi  to  P3,  aisi  and  bade  asimnth  from  Ps  to  Pi,  asi.  All  angular 
input  and  output  is  in  radians.  Compute: 

9i  =  tan~‘  [(1  -  f)  tan^,],  for  i  =  1, 2, 

AA  =  As  -  Ai,  C  =  (tfi  +  «a)/2,  =  (^s  - 

H  s  cos^  -  rin^  tfmt  ^  A^m  +  Psin^(AA/2), 
ds  2ain“^(I*/^),  17  =  28in*tfmCos*  A^m/(1  -  X), 

V’  =  28in2Atf«cQs3tfm/X,  X^U  +  V,  Y^U-V, 

T^df^d,  D  =  4T^,  ^  =  2cosi,  A  =  DE, 

C  =  T-{A-E)I2,  ni^X{A  +  CX), 

Y{B  +  EY),  ns  =  DXY,  (12) 

Sid  =  f{TX  -  r)/4,  S^d  =  /^(m  -  ns  +  n3)/64, 

5/0,  =  (r-5id  +  5jd) rind,  M  =  32T  -  (20r  -  A)Ar  -  (P  +  4)y, 

P  =  2y  -  P(4  -  X),  G  =  /r/2  +  /®M/64,  Q  =  -(PG  tan  AA)/4, 

AAi,  =  (AA+(J)/2, 

<1  3:qatn(-8inAtfmCOsAA|,,,coe^msinAA|,i), 
ts  =  qatn(cos  A^„  cos  AA;,  ,  rin  rin  AAJ,,), 

®i3  —  ti  +  ta  and  otsi  =  ti  “  ts* 

b.  Direct  Solution.  Input  variables  are  the  latitude  ^1  and  the  longitude  Ai  of 
Pi,  the  forward  asimuth  ais  from  Pi  to  P3,  and  the  distance  (5/a«)  from  Pi  to  P3.  The 
output  variables  are  the  latitude  ^  and  longitude  A3  of  P3,  and  the  backward  aiimuth 
031  from  P3  to  Pi.  All  angular  input  and  output  is  in  radians.  Compute: 

^1  =tan"^((l-/)tan^i|, 

M  =  cos^i8inoi3,  Cl  =  /M,  cs  =  /(I- Af^)/4, 

P  =  (1  —  C3)(l  —  C3  —  Cl  Ad),  P  =  C3(l  +  Cl  A#/2)/ Dy 
.¥  =  008^1008  013,  (Ti  =  qatii(¥,8intfi), 
d  =  (5/a,)/P,  tt  =  2(<ri  -  d),  =  1  -  2P  cos  u, 

V  =  co8(tt  +  d),  X  =  <^8indcoed(2y^  -  1), 

y  =  2PViy8ind,  A<r  =  d  +  X-r,  (13) 
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*  ■ 


K  —  [AP  +  [NcotAa  —  on^i  sin 
taa  ^  =  (sin  cof  A(r  +  JV  flin 

s  qatn(8in  A(r  sin  ai2 ,  cos  1  cos  Atr  -  sin  sin  Aff  coe  aia) , 

S  =  ei(l  -  ea]A^  -  CiCa  on A^ co8(20rt  -  A^), 

A3  »  Ai  +  Aq  — 

031  =qatn(-M,-(JircosA(r-sm9ismA<r)]. 

3.  Spherical  Sarth  Fommlaa.  The  inverse  solution  [Eqns.  (3),  (4)  and  (5)]  and  direct 
solution  [Equs.  (6),  (7)  and  (8)]  formulas  were  developed  in  Section  B  of  this  report.  They 
are  summarised  here  for  convenience. 

a.  Invene  Sohitkm.  Input  variables  are  the  latitude  ^1  and  longitude  Ai  of  Pi,  the 
latitude  ^  and  longitude  A3  of  P3,  ^  and  A3.  Output  variables  are  distance,  d,  forward 
asinmth  from  Pi  to  P3,  013,  and  back  asimuth  from  P3  to  Pi ,  031.  All  angular  input  and 
output  is  in  radians.  Compute: 

(i  =  arccosfsin^sin^i  +cos^coe^i  cos(A3- Ai)], 

013  sqatn(cos^ sin(A3-Ai), stuccos ^1- cos ^ sin ^icos( A3- Ai)),  and  (14) 
031  =qatn(-cos^isin(Ai  -  A3),8m^i  cos^ -cos^i8in^coe(Ai  -  A3)). 

b.  Direct  Sohitlom.  Input  variables  are  the  latitude  ^1  and  the  longitude  Ai  of  Pi, 
the  forward  asimuth  013  from  Pi  to  P3,  and  the  distance  d  Pi  to  P3.  The  output 
variables  are  the  latitude  ^  and  longitude  A3  of  P2,  and  the  backward  asimuth  031  from 
Pi  to  Pi .  All  angular  input  and  output  is  in  radians.  Compute: 

^  =  arcsin(co8d8in^i  +8indco8^i  coeoi3), 

A3  =  Ai +qatn(8indanoi3,cosdco8^i -sindsin^icoson),  and  (15) 
031  =  qatn(-  cos  ^1  sin  013,  sin  ^  sin  d  -  cos  ^1  cos  d  cos  013). 

e.  Given  Longitiide,  Find  Latitude.  Input  variables  are  the  latitude  ^  and 
longitude  Ai  of  Pi,  the  forward  asimuth  013  from  Pi  to  P3,  and  the  longitude  A  of  some 
point  P  on  the  great  circie  joining  Pi  and  P3.  The  output  variable  is  the  latitude  of 
P.  All  angular  input  and  output  is  in  radians.  Using  Equs.  (2),  label  the  angles  and  legs 
so  that  A  =  360*  -  ot3i ,  P  =  Qri3,  (7  =  A-Ai,a  =  90®-^i,5  =  90“-i  and  c  =  d. 


Snbfthriting  these  into  Bqne.  (2)  and  simpii^ring  the  diffeiencee  of  an^^  we  obtain 

cooaai  =  -coeaiacoe(A- Ai)  +  ainai3  8in(A- Ai)8in^, 
einaai  ein^  =  coeai3ein(A  -  Ai)  +  8inai3Coe(A  -  Aijsin^i,  and 
sin  031 001  ^  s  sm  ai3  coe  . 

Since  the  only  quantity  of  interest  is  and  since  latitude  lies  in  the  range  of  -90*  to  +90*, 
it  suffices  to  divide  the  second  equation  by  the  third  equation  to  obtain  an  expression  for 
tan  4’  We  obtiun 

-1  rcoe®i3ein(A  -  Ai) +8inoi3Cos(A  -  Ai)ain^i 

9  =  tan  - : - 7 - 

sin  013  cos 

Bbr  almost  all  ^iplkations,  Equ.  (16)  will  suffice.  It  is,  however,  possible  for  the  denom¬ 
inator  of  Equ.  (16)  to  become  lero  in  certain  circumstances.  We  now  face  a  situation  in 
which  qiecial  cases  must  be  examined. 

One  way  to  handle  the  dilemma  to  first  evaluate  the  denominator  and  the  numerator. 
If  the  value  of  the  denominator  is  sero,  then  ^  =  90*  if  the  numerator  is  positive  or 
4  s  ~oo*  if  the  numerator  is  negative.  If  the  denominator  is  not  lero,  then  use  Equ.  (16). 
A  second  way  to  handle  the  dilemma  is  to  compute 

4  -  qatn(co8  ais  8in( A  -  Ai )  +  sin  ai3  cos( A  -  Ai )  sin  ^i ,  sin  ai3  cos  ^i ] ,  (17) 

but  this  leads  to  another  problem  because  4  mny  now  be  in  the  range  of  —180*  to  +180* .  To 
correct  the  output  to  a  proper  latitude,  subtract  180*  if  ^  >  90*  or  add  180*  if  ^  <  -90*. 

The  question  of  which  method  to  use  is  a  matter  of  personal  choice.  Generally,  a 
simple  one  line  logic  fonction  can  handle  all  of  the  special  cases. 

d.  Find  the  Vertex.  The  vertex  of  a  great  circle  route  is  the  most  northerly  or/and 
the  most  southerly  point  on  the  great  circle.  Bowditch  [Ref.  1,  §1016,  Bqus.  B25-B30|  gives 
formulas  for  first,  computing  the  latitude  of  the  vertex,  and  second,  computing  the  longi¬ 
tude  of  the  vertex  as  a  fonction  of  the  latitude  of  the  vertex.  The  Bowditch  latitude  formula, 
^m  a  straight  forward  application  of  the  Law  of  Sines,  is  =  cos~^  (sin  ai3  cos  ),  where 
4^  is  the  latitude  of  the  vertex  and  4i  and  ais  have  the  same  meaning  as  in  the  previous 
section.  The  principle  angle  of  the  arccosine  lies  between  0*  and  180*  whereas  latitudes 
lie  between  -90*  and  +90*.  Bowditch  illustrates  only  cases  in  which  4*  i*  between  0*  and 
+90*  and,  unfortunately,  makes  no  mention  of  how  to  handle  cases  when  the  arccosine  is 
greater  than  90*. 


Once  A«  le  known,  ^  can  be  found  neing  one  of  the  methods  in  §e  such  ae  Eqn.  (17).  Since 
|A«  -  All  <  00*,  it  fdkws  that  ^  mnet  be  in  the  same  hemisphete  ae  if  is  in 

the  northern  (southern)  henuspboe,  then  most  be  the  northern  (southern)  vertex. 

e.  Ghrem  Lntitiide,  Pfaid  Losigitnde.  Input  variables  are  the  latitude  ^  and 
longitnde  Ai  of  Pi,  the  forward  asimnth  ai)  from  Pi  to  Pa  and  some  latitude  It  is 
convenient  but  not  necessary  to  know  the  latitude  of  the  vertex  of  the  great  circle  from 
Pi  to  Pa.  We  wish  to  find  the  lon9tnde{s)  on  the  great  circle  from  Pi  to  Pa  which  have 
latitude  Thoe  are  several  possible  outcomes  (1)  if  |^|  >  |^«|,  there  is  no  solution;  (2) 
if  ^  s  thm  is  one  solution  and  that  is  A  =  A«;  and  (3)  if  |^|  <  |^*|,  there  are  two 

longitudes  which  satisfy  the  given  conditions. 

Rewrite  Eqn.  (16)  in  the  form 

coeaiacos^8m(A- Ai)  +sinaia8in^icos^coa(A- Ai)  -sinaiacos^isin^, 
or 

5ain(A-Ai)  +  C’cos(A-Ai)  =  JT,  (20) 

where  S  s  oosaiacos^,  C  s  sinaiasin^i  coe^,  and  K  =  sinotiacos^  sin^.  Define  p>0 
and  tf  as  the  sohitions  to  the  equations 

Cspccmtf  and  S=:pmrf.  (21) 

Solving,  we  find  that  _ 

p  —  +  C*  and 

i»  =  qatn(5,C) 

Substituting  Equs.  (21)  into  Equ.  (20)  and  rearranging,  we  find 

coe(A  -  (Ai  +  ti)\  =  K/p. 

So  that 

A-(Ai  +7)  =  ±cos-‘  {K/p), 
or 

A  =  (Ai  + »»)  ±  cos”‘  {K/p) .  (22) 

Recall  that  a  knowledge  of  is  not  known  to  determine  the  number  of  solutions.  (1)  if 
\K/p\  >  1,  there  are  no  solutions;  (2)  if  \K/p\  =  1,  there  is  one  solution;  and  (3)  if 
\K/p\  <  1,  there  are  two  solutions. 


In  the  prerioiia  aectioii  it  wai  mentioned  that  there  waa  yet  another  WTCf  to  determine 
the  podtion  of  the  vertex,  lb  do  this,  for  fixed  define  /(A)  in  terms  of  ESqn.  (20)  as 

/(A)  =  58in(A  -  Ai)  +  Ccos(A  -  Ai)  -  K. 

Then  f(^),  the  derivative  of  /(A)  with  req)ect  to  A  is 

/(A)  *  5cos(A  -  Ai)  -  C'ain(A  -  Ai) 

Define  A*  so  that  /'(A*)  =  0.  Thns 

5cos(A*  -  Ai)  -  C8m(A*  -  Ai)  =  0, 


or 


rin(A*  -  At)  _  5 
cos(A*  —  AiJ  ”  O' 


A*  -  At  =qatn(5,(7), 

so  that 

A*  =  At  +  qatn(5,  C)  =  At  + 

Hence  A*  =  At  + is  an  extremnm  of  Eqn.  (20),  or  A«  s  A*.  Note  that 

/"(A*)  *  -58in(A*  -  At)  -  Ccoe(A*  -  At) 

=s  -5  sin  17  — (7 cos  7 

^  P  P 

S^  +  (P 


=  -/I  <  0. 

Thns  A«  s  A*  is  the  northeni  vertex. 

E.  Rhumb  Line  (Mercator)  Formolas. 

1.  Spheroid  Earth  Fornmias.  Bowditch  [Ref.  1,  Explanation  of  Thble  5,  Meridional 
Parts]  ^ves  the  formula  for  the  computation  of  meridional  parts  as 

M  =  a, b tan  ^45*  ^«^8b^  +  y  sb®^  +  ^ sb®^  +  . . 

=  Oebtan  ^45®  +  -  (ie«  ^esm^+  y  sb®^  +  y  ain®^H — ^ 


whoB  i/  b  the  munbcr  of  meridioiial  parti  between  the  equator  and  the  given  latitude 
4,  ii  the  equatorial  radiua  of  the  earth,  and  e  »  y/[2f  -  if  the  eccentricity  of  the 
earth,  where  /  s  1/298.26  ia  the  WGS  1972  earth  flattening  factor. 

Uiing  the  identity  [Bef.  7,  #766] 


it  ii  powble  to  rewrite  the  fonnnla  for  M  ai 


+  s)  (iris?) ' 


=  a«ln 


tan  (45-  + 1 


/l  +  edn^y 

^l-esin^  J 


If  we  denote  the  course  (forward  anmuth)  by  Oia,  then  [Bowditch,  VoL  II,  §1001  and 
§1013]  s  Aa  -  Ai  =  mtanoria  where  m  s  Af(^)  -  Af(h)t  that 


(Aj- AO/a, 


tan  46*  + 


/l  +  esm^y'*  /l  +  eein^iN 
^l-esin^y  ^l-«sin^iy 


This  equation  ia  also  given  by  Shufeld  and  Newcomer  [Hef.  8,  pp  81-8^.  It  ia  the  spheroid 
earth  rhumb  line  equation. 

e.  Larene  Sohitioii«  Input  variables  are  the  latitude  and  lon^tude  Ai  of  Pi,  the 
latitude  ^  and  longitude  A^  and  Pj,  ^  and  A^.  Output  variables  are  distance,  d,  forward 
asimnth  atg  from  P\  to  Pi,  and  back  aaimuth  a^i  from  Pi  io  P\.  All  angular  input  and 
output  is  in  radians.  In  Equ.  (23),  express  the  an^^es  in  radians,  then  to  obtain  an  in  the 
proper  quadrant,  rewrite  Equ.  (23)  in  the  form 


The  distance  is 


if  cosais  7^  0, 


<i=^  'cosflfia'  •  (2S) 

a«|  As  -  At  I  cos  ,  otherwise. 

In  the  WGS  1972  earth  model,  a«  =  3443.9174  nautical  miles.  Also,  an  =  an  +  t. 

b*  Direct  Sohitkm*  Input  variables  are  the  latitude  and  the  longitude  At  of  Pi, 
the  forward  asimuth  sis  from  P\  to  P),  and  the  distance  d  from  Pi  to  P^.  The  output 
variables  are  the  latitude  ^  and  longitude  As  of  Ps,  and  the  backward  asimuth  asi  from 
Ps  to  Pi.  All  angular  input  rnd  output  is  in  radians.  In  the  direct  rhumb  line  solution  for 
the  spheroid  earth,  if  cosai^  s  0,  then 

d/a« 


As  =  Ai  + 

<h  = 


cos^i’ 


and 


(26a) 


otherwise,  if  cos  oris  ^  0, 


As  —  Ai  +  tan  ois 


/l  +  eam^Y 
\^l-eain^y 
^  +  (d/a,)cosai3. 


In  tan  i 

/I  +  e8in^i  ^ 
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1  l-esin^i  ] 

and 


(266) 


Also,  otsi  »  ais  + 1. 


3.  Spherical  Earth  Fommlas.  The  spherical  earth  formulas  are  obtained  by  setting 
«  =  0  in  Equations  (24)  and  (26b),  and  replacing  a«  by  the  standard  approximation  of  60 
nautical  miles  per  degree  of  arc  on  the  earth's  surface. 

a.  Inverae  Solution.  Input  variables  are  the  latitude  and  longitude  Ai  of  Pi ,  the 
latitude  ^  and  longitude  As  and  Ps,  ^  and  As.  Output  variables  are  distance,  d,  forward 
asimuth  from  Pi  to  Ps,  ais,  and  back  asimuth  from  Ps  to  Pi,  asi-  All  angular  input  and 
output  is  in  radians. 

The  spherical  earth  rhumb  line  course  is  obtained  by  setting  e  =  0  in  Equ.  (24).  The 


result  is 


ois  =  qatn  ^As  -  Ai,btaa  -  In  tan 

Ibr  the  sphoical  earth  model,  the  distance  i  in  nautical  miles  is 


if  cos  ais  /  0, 


cos  ais 

1 60  |As  -  Aijcos^i,  otherwise. 


(27) 


(28) 


IS 


AjfOf  a)i  a  ai3  +  w. 

b.  DIraet  Sohttkm.  lapnt  variables  are  the  latitude  and  the  longitude  Ai  of  Pi, 
the  forward  asiniuth  Sis  frran  Pi  to  P3,  and  the  distance  d  from  Pi  to  P3.  The  output 
variables  are  the  latitude  ^  and  lon^tude  of  Pa,  and  the  backward  asimnth  a^i  from 
Pi  to  Pi.  All  angular  input  and  output  is  in  radians. 

If  cosaia  s  0,  then 


Aa  *  Ai  + 


[U0j[e0jami 


;  and 


(»*) 


otherwise,  if  cosaia  /  0, 


Aa  =  Ai  +  tanaia  j^lntan  j » 


(286) 


Also,  aai  »  aia  +  ir. 


IIL  NAVALGORi  Ncvigatkm  Algorithiii  Prograxn. 

A.  Xntrodiietioii.  NAVALGOR  Implementi  three  sete  of  compatationai  procedures:  (1) 
IKrect  Sointioa  Algorithms,  (2)  lavene  Solutm  Algorithms,  and  (3)  Rhomb  Line  Ap- 
pronrimations  to  Great  Circle  Routes. 

The  "Direct  Sobtion  Algorithms*  compute  the  latitude  and  lonptude  of  a  poeitioa 
and  the  backward  asimnth  from  Pj  to  Pi  given  the  latitude  and  longitude  of  a  position  Pi, 
the  forward  asimuth  from  Pi  to  Pi  and  the  distance  from  Pi  to  P^.  The  "frrverse  Solution 
Algorithms*  compute  the  distance  from  position  Pi  to  position  Pi,  the  forward  asimnth 
from  Pi  to  Pi,  and  the  backward  asimuth  from  Pi  to  Pi  given  the  latitude  and  longitude 
of  positions  Pi  and  Pi.  Fbr  compaiision,  the  direct  and  invene  computations  are  made 
using  four  procedures:  (1)  the  spheroid  earth  modri,  (2)  the  spherical  earth  model,  and 
rhumb  line  ^roodmations  using  both  (3)  the  spheroid  earth  and  (4)  the  spherical  earth 
models. 

The  spheroid  earth  model  should  be  considered  the  standard  of  comparison  for  all 
other  modds.  Computations  performed  mnng  an  IBM  soss  and  double  precision  FORTRAN 
to  compare  the  inverse  solution  algorithm  (Eq^  (12)]  to  the  agio  long  lines  (50  to  6000 
miles)  [Ref.  9|  demonstrated  that  the  tnagimnTn  average  error  of  the  inverse  solution  algo¬ 
rithm  was  -0.006  meters  with  a  standard  deviation  of  0.014  meters  (unpublished  result). 

The  spherical  earth  model  is  included  because  of  its  popularity  and  simplicity.  For 
many  procedures  the  spherical  earth  model  is  adequate,  and  it  certainly  is  more  compact 
for  use  in  a  small  computer  than  the  spheroid  earth  model 

The  "Rhumb  Line  Apprcodmations  to  Great  Circle  Routes*  procedure  may  be  used  to 
find  piecewise  constant  course  (rhumb  line)  approximations  to  great  circle  routes  for  any 
given  increment  in  longitude.  In  addition,  if  the  vertex  of  a  course  is,  for  example,  too  far 
to  the  north,  a  limiting  latitude  may  be  input  to  restrict  the  rhumb  line  approodmatiou  to 
go  no  further  north  than  that  limiting  latitude.  The  spheroid  earth  rhumb  line  equations 
are  used  in  this  section  of  the  program. 

In  the  sample  problems  and  in  the  program  listing,  the  east  minus  and  south  minus 
convention  is  used. 


B.  Sampk  Probkam  NAVALGOR*  The  ontpnt  shown  in  these  sample  problems 
was  computed  luing  Microsoft  BASICA/D  (double  precision)  on  an  IBM  PC. 

Muter  Mena.  The  master  menu  for  NAVALGOR  is  shown  below. 

MAVIGATIOM  ALGORITHM  DEMO 

1)  DIRECT  SOUmON 

2)  INVERSE  SOLUTION 

3)  RHUMB  LINE  AFPROmATIONS 

TO  GREAT  CIRCLE  ROinES 

4)  qurr 

Option  (1)  selects  the  *Direct  Solution  Algorithm*,  option  (2)  selects  the  ^Inverse  Solution 
Algorithm*,  option  (3)  selects  the  *Rhumb  Line  Approximations  to  Great  Circle  Routes* 
algorithm,  and  option  (4)  returns  the  user  to  the  operating  sjrstem. 

Problem  1.  Suppose  you  are  at  San  Francisco  (latitude  37* 47*  north  and  longitude 
122*26'  west),  that  your  initial  course  is  260*  and  that  you  travel  a  distance  of  4000  n.  mi. 
What  is  your  final  position?  Select  Option  1  ftom  the  master  menu. 

DIRECT  SOLUTION 

1st  UTmJDE  DD.MMSS  (-8)  7  37.47 
1st  LONGITUDE  DDO.MMSS  (-E)  7  122.35 
INITIAL  COURSE  DDD.MMSS  7  360 
DISTANCE  (n.  el.)  7  4000 

SPHEROID  EARTH  DIRECT  SOLUTION 
2Bd  UTITUDE  6*38.7' 

2aA  LONGITUDE  '172*08.8' 

BACK  AZIMUTH  51*40.7' 

SPHERICAL  EARTH  DIRECT  SOLUTION 
3Bd  UTITUDE  6*41.9' 

2iid  LONGITUDE  -172*00.7' 

BACK  AZIMUTH  51*35.9' 

RHUMB  LINE  SOLUnONS 

SPHERICAL  26*12.4'  -159*56.0' 

SPHEROID  26*12.4'  -160*18.4' 


PRESS  ANT  KET  TO  CONTINUE 


Probtem  2.  Suppoae  you  are  at  San  fVanciaco  (latitude  37*47*  north  and  lonjptnde 
122*25*  west)  and  that  your  destination  is  Sydney,  Australia  (latitude  33*51*  south  and 
longitude  161*13*  east).  How  to  do  you  travel,  what  is  your  initial  course,  and  what  is  the 
badcward  asimuth  from  Sydney  to  San  IVansisco?  Select  Option  2  from  the  master  menu. 

Z2VERSE  SQLimOH 

1st  UTnUDE  DD.MNSS  (-8)  7  37.47 
1st  LOMCnUDE  DDD.MN88  (-E)  7  122.26 
2ad  UTTITJDE  DD.NNSS  (-8)  7  -33.51 
2iid  LONGITUDE  DOO.MMSS  (-E)  7  -151.13 

SPHEROID  EARTH  INVERSE  SOLUTION 
DISTANCE  6443.52  a.ai. 

FORWARD  COURSE  240*29.3* 

BACK  COURSE  56*55.7* 

SPHERICAL  INVERSE  SOLUTION 

DISTANCE  6446.3  a.nl. 

FORWARD  COURSE  240*18.9* 

BACK  COURSE  56*46.9* 

RHUMB  LINE  SOLUTIONS:  SPHERE  SPHEROID 
DISTANCE  6464.49  6486.71  s.nl. 

COURSE  228*19.7*  228*29.7* 

PRESS  ANT  SET  TO  CONTINUE 

Problem  3.  Suppose  you  are  at  latitude  37*  north,  longitude  76*  west  and  that  your 
destination  is  latitude  45*  north,  longitude  1*  west.  Compute  the  initial  great  circle  course 
and  distance,  the  latitude  and  longitude  of  the  vertex,  and  a  rhumb  line  approximation  to 
the  course  traveUng  at  most  7*  degrees  of  longitude  on  each  leg  of  the  course.  In  addition, 
no  portion  of  the  course  is  to  be  more  northerly  than  47*.  Select  Option  3  horn  the  master 


RHUMB  Lm  APPROXIMAnOHS 


IHITIAL  UTTIUDR  X.MN88  (-S)  T  37 
iimAL  LoicrraDi  ddo.mnsh  (-»  t  re 
niAL  UmUDE  DD.MM88  (-S)  ?  46 
FIIAL  LOMCmiDl  D0D.MM88  (-K)  7  1 

OEAT  CIRCLI  SOLUnCi 

IIITIAL  C0UR81  SiVai.R' 

DISTAia  3307.8  a.al. 

RHUMB  LINI  (MERCATOR)  SQLUTIOM 
COURSE  81*58. a' 

DISTAMa  3436.0  a.al. 

VERTEX:  UTITUDE  48*10.8^ 

L0RGI7UDE  38*07.3^ 

PRESS  ART  KEY  TO  COHTIHUE 

In  this  fint  screen,  the  great  circle  solution  and  a  single  rhumb  line  solution  are  computed 
and  displayed.  In  addition,  the  location  of  the  nearert  vertex  is  displayed.  The  vertex 
may  or  may  not  lie  between  the  initial  and  final  positions.  If  the  vertex  is  not  between 
the  initial  and  final  positions  then  any  limiting  latitude  which  may  be  mput  on  the  next 

screen  prompt  will  be  ignored.  If  the  vertex  is  between  the  initial  and  final  positions,  a 

limiting  latitude  will  be  ignored  unless  it  is  lies  between  the  latitude  of  the  vertex  and  the 
latitude  of  the  position  which  is  closest  to  the  latitude  of  the  vertex. 

2nd  Sereenx 

RHUMB  LINE  APPROXIMATIONS 

INPUT  TEE  MAXIMUM  NUMBER  OP  DEGREES 
OF  LONGITUDE  ON  EACH  RHUMB  LINE  LEG  ?  7 

LIMITING  UTITUDE  DD.MMSS  (-S) 

(0  TO  OMIT)  7  47 
LIMITING  LONGITUDES: 

10*46.8'  46*28.8' 

PRESS  ANT  KET  TO  CONTINUE 

In  the  second  screen  a  prompt  appeared  for  the  maximum  longitude  between  each  rhumb 
line  leg — a  value  of  7*  was  input.  Also  input  was  a  limiting  latitude  of  47*,  which  lies 


between  the  vertex  latitude  of  48*19.8^  and  45*,  the  latitnde  of  the  final  poeition.  Limiting 
Imigitadee  of  10*45^  and  46*28^  are  computed  and  displayed.  1116  latitude  of  a  great 
drde  route  between  these  limitnig  lonptudes  will  be  more  northerly  than  47* ,  consequently 
any  rhomb  line  approximation  will  be  due  ^st  or  due  west  at  a  latitude  of  47*  between 
these  kmgitodes. 

Srd  Serecau 

RHUMB  LIRE  APPROXINAnOI  TO  GREAT  CIRCLE  COURSE 


GREAT 

GREAT 

RHUMB 

RHUMB 

GIRCLS 

CIRCLE 

LINE 

LINE 

UT 

LONG 

COURSE 

DIST 

COURSE 

DIST 

37*00.0' 

76*00.0' 

66 

0 

58 

0 

39*27.9' 

71*00.0' 

69 

278 

62 

279 

42*18.8' 

64*00.0' 

64 

639 

66 

641 

44*32.0' 

67*00.0' 

69 

971 

71 

974 

48*11.7' 

60*00.0' 

74 

1283 

78 

1287 

47*00.0' 

46*28.8' 

77 

1476 

90 

1480 

4roo.o' 

10*46.8' 

108 

2884 

105 

2900 

46*56,3' 

6*00.0' 

107 

3130 

108 

3148 

45*00.0' 

1*00.0' 

110 

3308 

108 

3326 

1)  NEW  APPROXIMATION  OR  2}  NEV  PROBLEM? 

Screen  3  displ^  the  rhumb  line  approximation  for  a  maximum  of  7*  of  longitude  between 
course  changes.  The  latitude  and  longitude  of  the  initial  and  final  positions  and  of  the 
positions  at  which  course  changes  are  made  are  given  in  columns  1  and  2,  respectively. 
Column  3  shows  the  great  circle  heading  at  each  poeition,  while  Column  4  shows  the 
cumulative  great  circle  distance  from  the  initial  poeition.  Similarly,  Column  5  shows  the 
rhumb  line  course  to  be  followed  between  each  pair  of  positions  and  Column  6  shows 
the  cumulative  rhumb  line  dirtance  from  the  initial  position.  Note  that  at  47*  north  and 
45*28.8^  west  (the  first  limiting  latitude  poeition)  the  rhumb  line  course  changes  to  due  east 
until  47*  north  and  10*45.8'  west  in  reached.  Also  note  that  the  rhumb  line  approximation, 
including  the  'detour'  at  the  limiting  latitude,  is  only  18  n.  mi.  longer  than  the  great  circle 
route.  The  third  screen  ends  with  two  options:  Option  1  returns  to  Screen  2  and  prompts 
for  new  inputs  while  Option  2  returns  to  the  master  menu. 


C.  Program  LktlBg— NAVALGOR. 


10  UN  ■lATIGinOH  ALG0IIIHN8"  I.H.  8HDD0I.  03-12-U.  BE?.  00-14-88  1600 
13  UN  "UTALGOB/BU* 

30  OEFOBL  A-T 

40  P4-Alir(l} ;P2*P4-»P4:PI-P2-»P2:TP-PI<^PI:IMI/180:EP-lE-33 
60  PL-l/2M.28:AM37813S/18S2:SC-flQA(n.«(2-n.)) 

60  OEFFni(Z)-X-MO«IHT(X/NO):BEN  Z  NOD  MO  PUNCnOIl. 

70  OEPFn.(Z)-Z-TP*IHT((Z«PI)/TP):BBi  LONOnOfiB  iDJUST  (-PI.PI) 

80  DEmB(Z)-IlfTa*N0«.8)/M0:IBI  ROONDSiO  FDMCTlOlt. 

90  0BPPW(Z)-Z<»PI»8GMaXAU(Z)>P2):UM  UnTOU  ADJUST  (-PI/2,PI/2) 

100  0EPFlfCa)-AIX(8qS(l’X*Z)/(Z-EP*a«<*)}}*PI*(Z<0i};UM  AICC08 
110  DEFFI8a)-ATX(Z/(8qB(l-Z*Z}-EP*(AB8CZ)-l))):UN  ABCSUf 
116  QOTO  2000 

120  IBI  oath  C-PI.PI)  PUHCnOM: 

130  A^ABI(T/(X-EP*(lii0#)))-PI*(X<0#)*(8GM(T)-a-0i))  :BEinU 
140  8B1  OATH  (O.TVQPl)  FUNCTION: 

160  A-ATS(T/a-EP*(X-Oi)))-PI*(X<Oi)+TP*(X>-Oi)*(T<Oi)  :BETUU 
170  BDf 

200  BEN  DIBBCT  SaumOH,  8PBEB0ID  EABTH.  ALL  ANGLES  MUST  U  IN  BA0IAN8. 

210  BOf  INPUT:  UHIUDE  01,  LONGITUDE  LI.  FORIABD  AZZNUTH  A1  AND 
220  BEN  DISTANCE  DI^<B/AE3  TO  A  POINT  P2.  NOTE:  DO  BAS  BADIAN  UNITS. 

230  BEN  OUTPUT:  UTITUDE  G2.  LONGITUDE  L2  AND  BACnABO  AZIXOTE  A2. 

240  8»4INai}:OCQS(Al):TA-ATN((l-PL)*TAN(01)) 

260  88-8IN(TA):C8-GQS(TA):N^80«C8:Cl-rL*M 
200  OPL*(1-M*N)/4:X1-G2)*(1-C2-C1*IO 
270  P«C2*(1<».6*C1«N)/D:N«G8*C0 
280  >N:X-88:G080B  130 :8G->A:8D>4»/D:IM* (80-80) 

200  N^1-2«P«C08(U)  :7M»8(U4>8D}  :X-C2«G2*8IN(8D)*C08(SD*(2*V*?-1)} 

300  T>2*P«V*N*8IN(8D):09-aD«X-Y:8L-8INC0q):CL-C08(Dq) 

310  I-8qB(N«lH-CN*CL-88*8L}*2)  :G2-ATN((88*a«N*8L)/X/(l-7L)) 

320  T— 8L«80:X«C8*<n.-B8*8L*C0:Ga8aB  1X:DJ-A 
330  B-Gl*(l-<:2}*Dq-Cl^*8L*C08(8G48G-Dq) 

340  L2-FNL(L1-^DJ-H}:T-N:Z>-(N*CL-88*8L):G080B  160:A2-A 
360  BETOU 
360  BEN 

400  BEN  INYEB8E  SOLUTION,  SPHEBOID  EABTH.  ALL  ANGLES  MUST  BE  IN  BADIANB. 
410  BEN  INPUT:  UTITUDE8  G1  A  G2.  AND  LONGITUDES  LI  A  L2. 

420  BEN  OUTPUT:  DISTANCE  0I>-(8/AE).  FOBVABD  AZDDTB  Al,  AND  BACXWABD 
430  BEN  AZIMUTH  A2.  NOTE:  DD  HAS  BADIAN  UNITS. 

440  TA-ATN( (l-FL) •TAN(G1} ) : TB-ATN( ( 1-PL) *TAN(G2} ) 

460  DN-L2-Ll:0T-(TB-TA)/2:TN*(TA<»TB)/2 

460  8B-8IN(DT)  :CH-C08(DT)  :8M-8IN(TM)  :C»<08(TM) 

470  H-CH*CH-8N*8N:L-8H*8H-»H*C8IN(DM/2))*2 

480  SlM^PNSCSqBCL)  ) :  U-2*8N«aM«GH*G8/  (  (l-L)  -  (EP*  (L-l)  )  ) 

400  7^*8H«U*01*CM/L :  X-U«7 :  T-U-7 :  T-8D/8IN(SD) 

SOO  I>4*T*T:B>2*C08(8D)  :A-D*E:0T-(A-E}/2:N1-Z«(A<K:«X) 

510  ►D+D:M>Y*(B*E*T)  :N>D*X*Y:D4-FL*rL*(Hl-M2+N3)/64 
620  D>PL«(T*X-Y)/4:D0-(T-D3-»04)*8IN(S0) 

630  M-32*T- (20*T-A) *X- (B*4) *Y : F»Y+Y-E*(4-X) 

540  G-FL*(T/2frL*ll/64}:q— (P*G*TAN(DN)}/4 
660  DV-(DM-»q}/2 :  8W«i8IN(DW) :  CXOSCDV) 


560  T4H*€«:X*<ai*8ff:G080B  130:T1-A 
570  I^-<a*Gi:X-8N«8W:GQ80B  iaO:T2-A 
580  M0-1P;Al-R0l(Tl<»T2}:i>rHN(Tl-n):ttTDIll 
500  UM 

000  in  oiucT  sQummi,  spiniciL  eixxb.  ill  AioLa  iibt  u  n  iadiais. 
610  m  hfut:  unTais  oi.  ldnoitude  li.  fouaid  azuoith  ai  aid 

680  UN  OUTAIRS  DO  TO  A  POUT  P2.  MOTE:  DO  18  H  AADIANS. 

630  UN  OOTFOT:  UTITUDE  02,  LONQnODE  18  HID  BAGXIAID  AZINOTH  A2. 

640  80-81X(A1):C»-C08ai):83-8IH(61}:C3-G08(61) 

660  8D-6IH(DD):CD<08(DD) 

060  T^8D*8O:Z-C3*GD-83*C0«6D:0080B  130:L2-nn.(Ll-A) 

670  7^-80*03 :X-flD*83-GD*CO*C3:G08UB  150:Aa-A 
660  G>n8(83<D«C3*80*ee):U7UU 

600  in 

700  UN  mVEUE  OOLOnOI.  8PHEIICAL  EAETH.  AU  AN6LE8  ND8T  U  IH  IADIAM8. 
710  UN  mPDT:  UTHUDU  61  A  02.  AHD  LONOITTIDEI  LI  A  L2. 

720  in  OUTPUT:  0I8XAHCE  DO  TO  A  POUT  P2.  (MOTE:  0  <-  DO  <-  PI  RADIAM8). 
730  in  fOEIAM  AZINOTH  Al.  AMD  BACZKAID  AZINDTI  A2. 

740  83-8IN(01)  :C3-Ca8(Gl)  :84-8IM(62)  :C4-C08(63) 

750  DL-L1-L2:8U-6IICDL}:CD-C08(DL} 

760  0>nC(83«844C3«C4«GU) 

770  >80*C4:XO*84-83*C4*CD:a08DB  150:A1-A 

780  T— 8D*<a:X-C4*83-84*G3*CD:G080B  150;A2-A:IZTTBM 

700  UN 

800  UN  UVEUE  lOLOnOM,  IHUNB  LIME  FOE  SPHEUCAL  A  SPHEIOTOAL  EAUI. 

810  UN  DtPUT:  UTITUDE8  01  A  02,  AHD  LONOITUDEI  LI  A  U. 

820  UN  OUTPUT:  0I8TAH(Z8  DA  A  DB,  AMD  AZIN0TH8  ZA  A  ZB. 

830  GAfTAM(P4-»01/2)  :GB-TAM(P4^2)  :DL-n(L(Ll-L3) : OI-An(DL/IO} 

840  OAFa-EP*(OAN))  :6B-GB-EP«(GB-0) 

850  OG»(02-01)/EO:T-OL:X«L06(GB)-LOG(GA):008DB  150:ZA^A 

860  CZ-008(ZA} :  0A-60*0E*C08(01} :  IP  AB8(a)> .  00017HEM  DAf60«0C/CZ 

870  E1«K*8IM(01) :  B>K*8IM(02}  :E>K/2 

880  >OL:X-LOO(OB/((ltE2)/(l-E2})‘ED)-LOO(GA/((ltEl)/(l>El))*ED} 

885  GQ80B  150:ZB"A 

800  CZ*G08(ZB):DB-60*DX«C08(01):1F  AU(a)>.0001THEH  OMO*DQ/CZ 
000  lETUU 

010  in 

1000  UN  DIIBCT  SOLUnOM,  UUNB  LINE  FOB  SPHEBICAL  A  SPHEBOIDAL  EABTl. 

1010  m  INPUT:  UniUDE  01.  LONGITUDE  LI.  FOUAID  AZIMUTH  Al  AMD 
1020  UN  DI8TANCE  DO  TO  A  POINT  P2.  NOTE:  DD  18  IN  BADIAN8. 

1030  in  OUTPUT:  UTITUDE  02  AND  LONGITUDE  U  A  L8 
1040  C6-C08(A1) 

1050  IF  06-OTHEN  LIpL1-»DD/C08(G1} :L8-U:G3-01:RETUU 
1060  02^1-»D0*C6 

1070  OA-TAN(P4<»01/2)  :GB-TAN(P4«02/2) :  TO-TANUl) 

1080  LIpL1-T6*(L06(GB}-L0O(GA)) 

1000  Bl-BC*8IN(01):E>BC*8IN(02):ED-BC/2 

1100  L8-U-T6*(L0G(GB/((1<»E2)/C1-B2})^}<LK(GA/((1-»B1)/(1-E1))*ED)) 

1110  lETUU 
1120  UN 

1200  UN  DECIMAL  TO  DDO  MN.F 

1210  •  ;IF  KOTHHI  V8-«-“;X— X 

1220  X-X^1/1200:Y-INT(X):7^7|+BIGHT|("  •♦8TRia) 


V * » * « * » *  *  "V.*  ■  1  r ■yy ."^. 


1230  M0O*a-T):y^Xna):2*-fTM(100O«Y) 

1240  V|-«|«ia])$(Z$.3.2)«V*-»UCSBn(Xt.l)«*'*:IElDU 
1210  IBI 

1200  UN  DDD.lOm  TO  DKINAL 

1270  n-0:7Qt  >1T0  LBtCVD.C^ilDKVI.Z.l):!?  O^.^TUH  XZ-Z 

1280  HECT  :Zr  IZ-OTBS  Z-7AL(Y|)  :SET0U 

1200  X-VAL(LEm(V|.ZZ)):nHl:ZF  Z<07HB  U— ra:Z— Z 

1300  V|-7|4"0000*  : >VALQIIM(V| .  n«l . 2)) : Z-?AL0lID|(7t .  IZ«3 . 2)} 

1310  Z-8V*((Z/60«Y}/60*Z}:UTDU 
1320  REN 

2000  CL8  :PU1IT  SPG(15):"NAViaiTZQX  ILOOftlTlN  OBA 

2010  PtINT  :PEIRT  :PSIHT 

2020  POINT  8PG(1S) ; "1)  DIRECT  SOUDTIQN 

2030  PUNT  8PC(1S)  :«2)  INTB8E  SQURIQN 

2040  PUNT  8PC(16}:*3}  RHOMB  LINE  APPROZINAnONB* 

2080  PUNT  8PC(1S):«  TO  GREAT  CIRCLE  ROOTEB 

2060  PUNT  :PUNT  8PC(15):"4}  QUH 

2070  G08DB  0010:C>?AL(C|} .ON  CG080B  3000.3500,4000.6500 

2080  GOTO  2000 

2000  REN 

3000  CL8  :PUMT  8PG(16);*DIRBCT  80L0T1QN":PUNT  :PUNT 
3010  PRINT  8PG(10);: PRINT  "18T  UTITDDE  DD.MOO  (-8)  *; 

3020  INPUT  7$:G08UB  1270:G1-Z*RD 

3030  PRINT  8PC(10}:: PRINT  ■18T  LaNGITOSB  DDD.MMU  (-E)  ■; 

3040  INPUT  V8:G08UB  1270:L1-X«U 

3060  PRINT  8PG(10):: PRINT  •INITIAL  COURSE  DDD.MMU  •; 

3080  INPUT  78:GQ8DB  1270:A1-Z*R0 

3070  PRINT  8PC(10)i«0I8TAMCS  (N.NI.)  "iiINPUT  D 

3080  01-0*RO/80:0>0/AE 

3000  G08DB  240:PRINT  :PRINT  8PC(8} ; •SPHEROID  EARTH  DIRECT  ULOnON 
3100  MI^100:G08DB  3210 

3110  01M)1:G080B  840:  PRINT  :  PRINT  8PC(8) ;  •8PBEUCAL  EARTH  DIRECT  SOLUnON 
3120  GOSOB  3210 

3130  GOSOB  1040:PRINT  :PRINT  8PC(8} ;  •RHUMB  LINE  ULOnOU*: 

3138  PRINT  8PC(4);^UT^:8PC(14):^L0NG 
3140  PRINT  8PC(12): •SPHERICAL  •; 

3160  Z-PN6(G2}:Z>iZ/RD: GOSOB  1210:V8<-V|:PRINT  7|;8PC(8} ; 

3180  Z-PNL(LR}:Z-Z/B0: GOSOB  1210:PRINT  Vf 
3170  PRINT  8PC(12): •SPHEROID  •;V8$;8PC(8}; 

3180  Z-PNL(LS};I-Z/R0:G08DB  1210:PRINT  VI 
3100  GOSOB  8000: GOTO  2000 
3200  REN 

3210  PRINT  8PC(13},^2ND  UTITODE  •; :I-03/RD: GOSOB  1210: PRINT  V| 

3220  PRINT  8PG(12):^2ND  LONGITUDE  •; :Z-L2/RD: GOSOB  1210:PRINT  V| 

3230  PRINT  8PC(12}:*BAa  AZIMOTH  •;  :X-A2/RD:G080B  1210 -.PRINT  V|:RETUU 
3240  REM 

3800  CL8  :PRINT  8PC(18) i^INTERU  SOLOTIQM*: PRINT  ;PRINT 
3610  PRINT  SPCCIO);: PRINT  •IST  UTITODE  DD.MISS  (-8)  •; 

3820  INPUT  T|:G08UB  1270:61-X 

3830  PRINT  SPCCIO):: PRINT  •IST  LONGITUDE  DDD.MMU  (-E)  ■; 

3640  INPUT  V|:G08UB  1270:L1-X 

3660  PRINT  SPCCIO);: PRINT  *280  UTITODE  DD. 10(88  C-8)  •; 

3880  INPUT  V|:G08UB  1270:G>X 


24 


3870  PBZIT  8PCao)::PftIHT  >3110  LOMnODB  000.1008  (-B)  ■; 

3880  UFOT  78:G080B  1370:I.>X 

3800  ai-01*ID:a>a2*IO:Ll-U*80:L2-U*80 

3000  QQ8QB  440:PEI]IT  :PUHT  8PG(8) ; >8PHEI0Z0  BAITI  ZMVEI8E  SOUmOV 
3610  ll>100:n-l:G08IB  3700 

3620  00-01:00800  740:P1I1IT  tPftZXT  8PC(8) ; >8PBmCAL  DIVBtfS  SQLOTIQX 
3630  lX-a:G08DB  3700:00800  830 

3640  PUHT  :ntlHT  81>G(8) ; >10010  LZHB  OOLOTIONO:  SPBEIB  8PHEI0I0 
3680  P8ZBT  8PG(ia) ; >01878808  >:F110(Di}:8PC(8): 788(08) ;>  8.81. 

3660  P8Z8T  870(12} ; >000888  >; :X-ZA/80: 00800  1210:78187  V8:870(6); 

3670  Z-Z8/80:00808  1210:78187  Vl 
3680  00800  6000:0070  2000 
3600  888 

3700  78187  870(12) ; >01871808  >; 

3710  U  iZ-ina  P8IB7  FmCDD*W:‘  N.Nl. 

3720  17  1X-27BD  78187  788(60*00/80) ; >  8.81. 

3730  78187  870(12) ; >7088180  000888  >; :Z«11/BD:0080B  1210:78187  78 
3740  78187  870(12); >8108  CQ0888  >; :Z-12/BD:0Q80B  1210:78187  78:887088 
3780  BBl 

4000  CL8  :78187  870(16); >88088  L18B  1P780Z18171088>: 78187  : 78187 
4010  78187  870(10) ; :78ia7  >181711L  L1717DD8  00.8888  (-8)  >; 

4020  18707  71:00808  1270: 01-1*80: 06-01 

4030  78187  870(10) ; :78187  >1817111  L0801700I  000.18188  (-8)  >; 

4040  18707  78:00808  1270:L1-Z*B0:L^L1 

4050  78IH7  870(10) ;  :78IH7  >n81L  UTllOOE  00.8888  (-8)  >; 

4060  18707  78:00808  1270 :a>Z*80: 07-02 

4070  78187  870(10);  :781B7  >7I81L  LOMOllOlS  D00.M888  (-B)  >; 

4080  18707  78:00808  1270: 12-1*80 :L7-L2 

4000  00808  740:80-10:78187  :  78187  870(10) ;>QUE17  dBCLE  80L07I0H 
4100  78187  870(16):  >1817111  000888  >;  :1>11: 1-11/10:00808  1210:78187  78 
4110  78IB7  870(18):  >01871801  >; 788(60*00/80) ;  >  N.NI. 

4120  0Q80B  830:78187  : 78187  870(12); >88080  L18B  (88801708}  80LU7108 
4130  78187  870(15): >000888  >::X-ZB/B0:Q080B  1210:78187  78 
4140  78187  870(16): >01871808  >:r88(DB);>  N.81. 

4160  G080B  4610 

4160  7BI87  : 78187  870(12); *788781:  UlITOOB  >; :Z-07/1D:G080B  1210:78187  78 
4170  78187  870(20) ;>L0MGI7O0B  >; :Z-L7/B0:GQ80B  1210:78187  78 
4180  GOSQB  6000 

4190  CL8  :78I87  870(15);  >88088  L188  1P7BOZD(1TI088>: 78187 
4200  78187  :78in7  870(10); >18707  TBB  81X1X08  ND8BQI  OF  DECREES 
4210  78187  870(10):*  OF  LONCITOOE  ON  EICH  RB08B  LINE  LEG  >; 

4220  18707  78:GQ80B  i270:I8*iX*BD 

4230  G8-0:B8>«:IF  I8-0THER  I8-T7 

4240  DL-P8L(L7-L6):1F  OL-OIBEN  F-1 

4280  07^F8L(L7-L8):IF  DLOOTHEN  F^D7/0L 

4260  LX-0:IF  F<-001  F>-17HEN  11-1:0070  4370 

4270  78187  :  7B187  870(10) ;  >LINI7I8G  UniODE  00.8888  (>8) 

4280  71187  870(10);*  (0  TO  0817)  *;:I8707  7|:G080B  1270:L>X*RD 
4200  L8-0:  L>0 : 61-G6 :  G3-G1 :  L1-L6 :  Ul-Ll :  7H-LL :  Ga-C7 :  L2-L7 
4300  IF  LL-07HDI  11-1:0070  4370 

4310  91-188(77-06) :  q>lB8(77-07)  :l-qi :  IF  q2<l  TBEH  A-q2 
4320  1-1*808(77)  :>7V-LL:  IF  1-0  THEN  0-6:0070  4340 
4330  >6/1 
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4340  IF  OC-O  01  a>l  im  LXf1:P1ZMT  8PC(10):"CA1IT  UIB.  -  ZairaiED.«:GOTO  4370 
4310  amnff  4810:P1IHT  8PC(10):"LI1IITIH0  LQBeiTDOES: 

4300  3^L8/ID:Q0iaB  iaiO:PllHT  8PC(ia):V|:8PG(10):X-LO/BD:008aB  iaiO:PlIFr  V| 
4370  60803  6000:60803  6230 

4380  0Mn(I.8-Ll):0^nn.(LB-Ll):ZF  138(00) <438 (D8) THEM  T-L0:L>L8:L8-T 

4300  lQFl:LX-L6:PX-66:Zr  I>1THE8  4460 

4400  L>L8:FT-LL: 60803  4010 

4410  U«L8:61-LL:I.2-L0:63-LL:608DB  740:60803  830 

4430  60803  8110:68-68«DD:I8-I8<»D3 

4430  IJ-L0:P>LL:L1-LZ:01-R 

4440  L2-L7:63-«7 

4460  NX-0: LT^L7:PT^: 60803  4010 

4460  lEM 

4470  PUHT  :P1INT  8PG(10}:*1}  MEN  APPB0ZIM4TIQN  01  2)  HBT  PIOBLZMT 
4480  60803  8010:G-VAL(C8):a>  C608DB  4100.2000 
4400  6070  4460 
4800  BBI 

4600  8EH  CONFOTE  VEITEX 

4610  81rtZX(lZ):a-GQ8(lZ}:8P-8IN(61}:CP-G08(61) 

4620  17  8A-0nn  67^72:  LV^L8:1ETDRN 

4630  L7^Ll-lTN(l/8P/TiNClZ)):U-Ll:LV^7SL(LT):LN>LV:6080B  4720 
4640  6?^:IEX0IM 
4660  lai 

4700  ROi  FIND  UTIinOE  6IVE6  LONOITDDE 

4710  81-8IX(lZ}:a-GQ8(lZ):8P-8IN(Pl):CP-C08(PA) 

4730  DL-U-UI:1MP«G08(SL)«81^8I1I(DL}*G1:Z-CP*81:6080B  130: Pi-1 
4730  P»4i6(P8):BEiaiM 
4740  BBI 

4800  BBI  FIHD  LOHGITODE  6IVBI  UTITDDS 

4810  IX-0:ZF  1B8(PB)->AB8(67)TBEII  KX-1:BET0IN 

4820  8MIB(lZ}:a-C08(lZ):8P-8IN(6a):a>-C08(63) 

4830  I><P*81*TlM(PH):C-8P481:8-Gl:B0-8qi(8*8K:*C}:Y-8:Z-C:608DB  130:NI>1 
4840  CC-fiC(Z/B0):DD-Lll-ini:L8-FHL(DD-0C):L»>Fin.(]B»GC}:BETDBN 
4860  BEN 

4000  BEN  BHDNB  IPPBQZZNinaNS  BUBBOnTIliS.  mPOT:  LX.  PZ.  LT  8  PT. 

4010  OL-FHLCLT-LZ) :  ZC-AB8(IX) :  16-138 (DL)/IC 
4020  IF  lC/BO<lTHDr  IIP-DITCl6+.5):MO-.l 
4030  IF  IC/BO>-lTHDi  NP-INT(16}  :M^1 
4040  IF  NP<1TIEN  IODL:L>LT:6QTa  4070 
4060  IC»86N(DL}*IC 

4060  L2-nm(  (U-»(DL-  CRP<  1)  *16}  /2}/BD-» .  6«86HCDL) )  *BD 

4070  L1«LZ:61«PX: 111-12:60803  4720:63-78 

4080  60803  740:60803  830:60803  6110:G8-Q8-»DD:B8-R8’»DB 

4000  IF  NP<1THBI  6060 

6000  NC-1:IF  NP-ITHEN  6040 

6010  L1-L2 : 61-62  :L2-PNL(L1’»IC}:L1HL2:6080B  47aO:6>Pl 
5020  60803  740:60803  830:60803  6110:68-G8-»D0:B8«R8*1» 

6030  )iC-NCfl:IF  NC<NPTHEN  6010 
6040  L1-L2: 61-62 :L2-LT:62-PT 

6060  60803  740:60803  830:60803  6110:68i<8-»DD:B8-B8«D3:IF  MX-lTiBf  BETDBN 
5060  U-L7:01-PT:liO-1P;ll-Fiai(12«PI}:G080B  6110;BETDBN 
6070  BBI 

6100  BBI  PBIHT  LINE  OF  OOTPOT 


BllO  l»il:Ml/ID:QC»OB  1210:niHT  Vl:8PC(2): 

S120  Ml/lD:a0S0B  1210:PiZ>T  T|; 

Biao  >nil(Al/IO):a(»IB  S200:PBIMT  y|:8PC(4}i 
6140  X-nii(eO*OB/IO)  :GOSQB  6200:PftZHT  Vt:8PC(4}: 

6160  X-Pn(ZB/IO}:Q08DB  6200:nZ>T  T$:8PC(3}: 

6100  Z-m(ll):G08aB  6200:PtIHT  ?| 

6170  n.-n.<»l:lF  ML-Wnn  OOSUB  OOOOzOQOQB  6220 
6180  mUM 
6100  UM 

6200  "•»8n8(Z):V|-ftZflnt0n.7):ISRnUt 
6210  BBC 

6220  08  :PtZ8T  8PC(10};*ina  LSB  iPPtOZllUnmi  TO  GBEiT  CJBOZ  CODISE*: 
PBXIT 

6230  PBBIT  *  6EE&T  GIEiT  IBDNB  RHOMB 

6240  PEIHT  ■  CnOB  CXBCLE  LIHE  LdB 

6280  PIXHT  ■  UT  LDHO  CQ0R8E  OUT  COUBSB  DI8T 

6200  IfL-6:lETDIH 

6270  IBM 

6800  08  :Bn) 

6810  RBI 

0000  PlZn  :PR1HT  8PG(10):*PRE88  AMT  OT  TO  CQMTXMDZ 
0010  FOE  ZI-ITO  g:C|-ZMIET8  :Mm 
0020  C^mm  :IF  C^««T1BI  0020 


D.  Program  Aimotatkm — ^NAVALGOR. 


LUM(f) 


140 

160 

240-350 

440-580 

640-680 

740-780 

830-900 


1040-1100 


1200-1240 

1270-1310 

2000-2060 

2070 

3000-3230 

3080 


UsagB 

PI4  =  ir/4.  PB  =  ir/2.  PI  =  ir.  TP  =  2ir.  HD  =  ir/180  is  the  degree-to-radian 
• _ 

oonvenioQ. 

FLls  the  earth’s  flatteaing  factor.  A£  is  the  earth’s  equatorial  radius.  (Both 
of  these  constants  aze  ftom  the  WGS-72  earth  model.)  ECC  is  the  earth’s 
eccentricity. 

FNM  is  the  X  mod  MO  fnnction. 

FNL  adjusts  longitude  to  lie  betwemi  -x  and  x. 

FNR  rounds  z  to  the  nearest  MOth. 

FNG  adjusts  latitudes  to  lie  between  -x/2  and  x/2. 

FNACS  is  the  arccoeine  fnnction.  In  the  Sharp  P^ISOOA,  fnnction  calls  may 
be  replaced  inline  by  the  ACS  function. 

FNASN  is  the  arcsine  fnnction.  Li  the  Sharp  PC-1500A,  fnnction  calls  may  be 
replaced  inline  by  the  ASN  function. 

FNATN2  is  the  qatn  function  whkh  returns  a  principle  value  between  -x  and 

X. 

FNATNP  is  a  qatn  fnnction  which  returns  a  principle  value  between  0  and  2x. 

GOTO  2000  to  commence  execution. 

Computation  of  the  spheroid  earth  direct  solution.  G1  s  di*  ^1  =  ^i>  ^1  = 
aia  and  distance  DD  are  input.  G2  L2  =  Aa  and  A2  ^  aai  are  output. 

C<nnputation  of  the  spheroid  earth  inverse  solution.  Gl  ss  di,  Ll  s  Ai,  G2  s 
^  and  L2  s  Aa  are  input.  A1  =  aiat  A2  =  oai  and  distance  DD  are  output. 

Computation  of  the  spherical  earth  direct  solution.  G1  s  Ll  s  Ai,  Al  » 
aia  and  distance  DD  are  input.  G2  =  L2  =  Aa  and  A2  =  aai  are  output. 

Computation  of  the  spherical  earth  inverse  solution.  G1  =  di>  LI  =  Ai,  G2  = 
^  and  L2  =  Aa  are  input.  A1  =  oia,  A2  —  oai  and  distance  DD  are  output. 

Computation  of  the  rhumb  line  inverse  solution.  Gl  =  dii  Ll  =  Ai,  G2  =  ^ 
and  L2  =  Aa  are  input.  Sphoical  earth  distance  DA  and  asimuth  ZA  as  well 
as  spheroid  earth  distance  DB  and  asimuth  ZB  are  output. 

Computation  of  the  rhumb  line  direct  solution.  Gl  =  Ll  =  A],  A1  = 
oia  and  distance  DD  are  input.  G2  =  spherical  earth  longitude  LR  and 
spheroid  earth  bngitude  LS  are  output. 

Convert  decimal  degrees  to  degrees,  minutes  and  tenths  of  minute  format. 

Convert  packed  degrees,  minutes  and  seconds  format  (DDDJ4MSS)  to  decimal 
degrees. 

lldaster  menu  option  display. 

Transfer  to  seized  menu  item. 

Direct  solution  option. 

D  is  the  distance  (input)  in  nautical  miles,  D1  is  the  spherical  earth  distance 
in  radians  and  DD  is  the  distance  in  earth  equatorial  radius  units. 

Compute  spheroid  earth  direct  solution.  Print  heading. 

Print  spheroid  earth  direct  solution.. 
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3110  Compute  spherical  earth  direct  solution.  Print  heading. 

3120  Print  spherical  earth  direct  solution. 

3130  Compute  rhumb  line  solution.  Print  heading. 

3140^160  Print  q)herical  earth  rhumb  line  latitude  and  longitude. 

3170-3180  Print  ^)herokl  earth  rhumb  line  latitude  and  longitude. 

3100  Continue  prompt.  Go  to  master  menu. 

3210-3230  Output  subroutine  for  option  1. 

3500-3740  bivm  solution  option. 

3600  Compute  spheroid  earth  inverse  sdution.  Print  heading. 

3610  Print  spheroid  earth  inverse  solution. 

3620  Compute  spherical  earth  inverse  solution.  Print  heading. 

3630  Print  spherical  earth  inverse  solution. 

3640  Heading  for  rhumb  line  solutions. 

3650  Print  q)herical  and  spheroid  earth  rhumb  line  distances. 

3660-3670  Print  ^)herical  and  q)heroid  earth  rhumb  line  courses. 

3680  Continne  prompt.  Go  to  master  menu. 

3700-3740  Output  subroutine  for  option  2. 

4000-5260  Ohnmh  line  approximation  to  great  drde  option. 

4000-4080  Input  prompts. 

4000-4110  Compute  great  circle  solution  and  print  initial  course  and  distance. 

4120-4140  Compute  rhumb  line  solution  and  print  course  and  distance. 

4150-4170  Compute  and  print  the  latitude  and  longitude  of  the  vertex  of  the  great  circle 
route. 

4180  Display  continue  prompt. 

4190-4220  Prompt  for  the  longitude  increment  of  the  rhumb  line  approxiination. 

4230  GS  and  RS  are  the  cumnlative  great  circle  and  rhumb  line  distances  traveled  on 
each  leg,  respectively.  Initialise  them  to  lero.  IN  is  the  longitude  increment 
in  radians.  U  input  as  sero  (no  increments  requested),  it  is  set  to  2%. 

4240-4260  Determine  if  the  vertex  lies  ^tween  the  origin  and  destination.  If  not,  go  to 
4370.  Otherwise  proceed. 

4270-4280  If  the  vertex  is  on  the  great  circle  route,  prompt  for  a  limiting  latitude  LL. 

4290-4340  Set  limiting  longitudes  L8  and  L9  equal  to  sero.  Determine  if  the  limiting 
latitude  cuts  the  great  circle  course.  If  it  does  not,  inform  the  user  that  the 
limiting  latitude  is  ignored.  L%  is  sero  if  the  limiting  latitude  is  to  be  used, 
otherwise  it  is  one. 

4350-4360  If  the  Umiting  latitude  cuts  the  great  circle  course,  compute  and  display  the 
longitudes  L8  and  L9  at  whidi  the  limiting  latitude  cuts  the  great  ciiole. 

4370  Display  continue  prompt,  then  print  heading. 

4380  Determine  which  limiting  longitude  is  closest  to  the  initial  position.  L8  becomes 
the  closest,  L9  the  farthest. 

4390  LX  and  PX  are  the  initial  great  circle  longitude  and  latitude.  If  the  limiting 
latitude  is  not  to  be  used,  go  to  4450. 
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4430-4440 


4470-4400 

4610-4640 

4710-4730 

4810-4840 


4810-5060 


4910-4050 


4060-4080 

4000 

5000 

5010-5030 

5040-5050 

5060 

5110-5200 

5220-5260 

5500 

6000-6030 


LY  and  PY  are  the  longitude  and  latitude  of  the  final  point  on  the  first  segment 
of  the  great  circle  course.  GOSUB  4810  to  compute  the  rhumb  line  approxima¬ 
tion  up  to  the  first  limiting  lonptude. 

Compute  the  great  circle  and  rhumb  line  distances  from  the  limiting  latitude 
at  the  first  lifting  longitude  to  the  second  limiting  longitude. 

Print  one  line  of  output  (for  the  limiting  longitudes).  Update  the  distance 
counters. 

Reset  LX,  PX,  LY  and  PY  for  the  section  of  the  great  circle  and  rhumb  line 
following  the  limiting  latitude  section  of  the  course. 

GOSUB  4810  to  compute  the  final  section  of  a  course  following  a  limiting  latitude 
leg  or  compute  the  entire  course  for  cases  in  which  the  limiting  latitude  is  not 
us^. 

Prompt  for  either  a  rework  of  the  rhumb  line  approximation  with  new  param¬ 
eters  or  for  a  new  problem. 

Compute  the  latitude  GV  an  longitude  LV  of  a  great  circle  vertex. 

For  a  given  longitude,  compute  the  corresponding  latitude  of  on  a  great  circle. 

Fbr  a  given  latitude,  compute  the  corresponding  longitudes  L8  and  LO  on  a 
given  gr  .it  circle.  is  set  to  one  if  there  is  no  solution,  otherwise  K%  is 
aero. 

This  is  the  subroutine  which  selects  the  longitude  endpoints  of  the  rhumb  line 
approximation  to  a  great  circle  between  a  starting  position  with  latitude  PX 
and  longitude  LX  and  a  final  position  with  latitude  PY  and  longitude  LY.  It 
then  coii^)utes  and  prints  the  heading  and  cumulative  distance  along  each  leg. 

Using  the  mcrement  sise,  IN  and  IC,  compute  the  number  of  interior  points, 
NP,  on  the  rhumb  line  approximation  exclusive  of  the  initial  and  final  points 
(the  number  of  interior  legs  is  NP  -  1.  Redefine  IC  if  it  is  longer  than  the 
difference  of  the  initial  and  final  longitudes. 

Find  the  heading  and  distance  on  the  initial  leg.  Print  out  first  line  of  output. 

If  there  are  no  internal  legs,  go  to  5060. 

Initialise  the  counter  1C.  K  there  b  only  one  leg  remaining,  go  to  5040. 

Loop  to  compute  print  each  internal  leg. 

Compute  the  headings  and  distances  for  the  final  leg.  If  is  one,  return  to 
compute  the  limiting  latitude  leg. 

Print  the  final  summary  line  of  output. 

Subroutine  to  print  one  line  of  output.  NL  counts  the  number  of  lines  of  output, 
if  20  lines  are  output,  a  new  screen  is  started. 

Subroutine  to  print  screen  heading. 

Return  to  BASIC. 

“PRESS  ANY  KEY  TO  CONTINUE’  subroutme. 


IV.  NAVEPHMx  Ahnanac  and  Ephameria  Program. 

A.  Introduction.  The  basis  for  NAVEPHM  are  the  equations  of  VanFlandem  and  PuUcki- 
nen  (VFP)  (Ref.  10|.  These  equations  can  be  used  to  compute  the  mean  heliocentric 
positions  of  the  sun  and  planets  and  the  mean  geocentric  position  of  the  moon  for  the 
mean  ecliptic  and  equinox  of  date.  The  anthon  claim  their  formulas  to  be  of  low-precision 
(!')  and  valid  for  any  epoch  within  300  years  of  the  present.  When  corrected  for  nutation 
and  aberration  the  accuracy  of  their  formulas,  at  least  for  the  sun,  moon  and  navigational 
planets,  appears  to  be  much  better,  i.e.  0.2^,  at  least. 

The  formulas  of  VFP  are  used  to  compute  heliocentric  spherical  ecliptic  coordinates 
for  any  specified  ephemeris  time,  ET.  These  coordinates  are  longitude  A,  latitude  and 
radius  vector  r.  These  coordinates  must  be  converted  to  rectangular  coordinates  x,  y,  and 
z  using  the  standard  transformation 

X  =  rcos/9cosA, 

y  =  rcos/JsinA,  and  (30) 

z  = 

To  obtain  the  geocentric  x,  y,  and  z  coordinates  of  the  planets,  subtract  the  x,  y,  and  z 
coordinates  of  the  sun  firom  the  x,  y,  and  z  coordinates  of  the  planets,  respectively. 

The  Julian  day  number,  required  in  many  calculations,  is  obtained  using  the  equation 
on  page  B2  of  Ref.  2.  The  ephemeris  time  is  obtained  from  the  universal  time,  UT,  from 
the  equation  ET  =  UT  +  AT.  The  factor  AT  is  obtained  by  astronomical  observation 
only.  The  formula  used  here,  AT  =  81.94T  —  15,  was  obtained  by  least  squares  from  the 
1980.5  through  1985.5  values  given  on  pages  B5  and  K9  of  The  Astronomieal  Almanac 
1985  [Ref.  11],  where  T  is  the  number  of  Julian  centuries  elapsed  from  1900  January  0, 
12^  ET.  Although  a  plot  of  AT  as  a  function  of  time  is  linear  for  1980.5  through  1985.5, 
this  should  be  checked  with  each  new  edition  of  The  Astronomical  Almanac.  An  accurate 
value  of  AT  affects  only  the  computation  of  the  moon’s  position.  Errors  of  as  much  as 
9.6*  in  AT  will  affect  the  computation  the  moon’s  position  by  at  most  0.1',  consequently 
it  will  not  be  necessary  to  change  FMPL  distribution  tapes  yearly. 

The  heliocentric  in-plane  velocity  components,  ^  and  yu .  of  the  planets,  required  for 
the  aberration  correction,  can  be  computed  from  the  formulas  [Ref.  12,  pg.  85] : 

Xw  =  -oeEsa.  E,  and 
Hu  =  aEy/ 1  -  cos  E, 


(31) 


wliere  E  is  the  eccentric  anonuly,  ^  is  the  time  derivative  of  a  is  the  semimajor  axis  and 
e  is  the  eccmtricity.  Unfortnnately,  to  find  E  one  must  solve  Kepleris  equation  iteratively, 
which  is  a  abw  process.  With  an  error  of  no  more  than  2%  for  the  navigational  planets, 
one  can  use  the  mean  anomaly  Af  in  place  of  £  in  Bqus.  (31).  The  in-plane  velocity 
components  are  thus  ^proodmated  by 

Cm  =  -oe^smM,  and 
%  =  ai?>/l-e*cosM. 

fVom  Eqn.  (3.76)  of  Ref.  12,  it  can  be  shown  that 

where  is  the  reduced  mass  and  k  is  the  gravitational  constant.  The  heliocentric  ecliptic 
velocity  components,  2,  y  and  i,  can  then  be  obtrined  by  the  transformation 

'x1  fcosn  -anO  Ol  fl  0  0  1  Tcosh;  -sino;  0]  Tzm 

y  sinfl  coed  0  0  cost  -sin*  sinw  cosw  0  Vw  , 

i  0  0  10  sin*  cost  0  0  10 

where  0  is  the  longitude  of  the  ascending  node, « is  the  orbital  inclination  and  u  is  the 
argument  of  perihriion.  The  equations  for  computing  a,  e,  M,  u,  *  and  d  were  obtained 
from  Escobal  (Ref.  13,  pp.  8-9].  The  geocentric  x,  y  and  z  coordinates  of  the  planets  are 
obtained  by  subtracting  the  x,  y  and  i  coordinates  of  the  sun  from  the  heliocentric  x,  y 
and  i  coordinates  of  the  planets.  Since  the  aberration  of  the  moon  is  negligible,  its  velocity 
components  are  not  computed. 

Ouce  the  mean  geocentric  rectangular  positions  (xm,  ym  and  z^)  and  velocities  (xm, 
ym  and  im)  have  been  obtsJned,  the  longitude  and  latitude  for  the  mean  ecliptic  and 
equinox  of  date  can  be  obtrined  by  inverting  Equs.  (30),  so  that 


A„  =  qatn(ym,x,»),  and 


Next  the  nutation  of  longitude  and  obliquity  Ac  are  computed  [Ref.  4,  §2C|.  The 
geocentric  rectangular  positions  (xt,  yt  and  zt)  for  the  true  ecliptic  and  equinox  of  date 
can  be  obtained  from  the  transformation 

’xt'  CcosA^  -sinA^  0]  ["xm 
yi  =  sinA^  cosA^^  0  ym  . 

[  0  0  ij  [zm 
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Tlie  same  tiansfonnation  is  osed  for  the  velocity  components.  The  longitnde  and  latitude 
for  the  true  ecliptic  and  equinox  of  date  can  be  obtained  by  substituting  the  true  positions 
into  Equs.  (33). 

The  aberratioa  or  light-time  correction  converts  true  positions  into  ^parent  positions. 
The  x-coordinate  transformation  is 


rxt 

, 

c 

where  r  is  the  geocentric  distance  and  e  is  the  velocity  of  light.  Similar  transformations 
are  used  for  the  y-  and  s-coordinates.  The  apparent  lon^tude  and  latitude  for  the  true 
ecliptic  and  equinox  of  date  can  be  obtained  by  substituting  the  apparent  positions  into 
Equs.  (33). 

The  apparent  geocentric  rectangular  positions  (xo^,  and  Sa«)  for  the  true  equator 
and  equinox  of  date  are  obtained  from  the  transformation 


1  0 
0  cosAe 
0  sinAe 


0 

■sinAe 
cos  As 


The  apparent  ri^t  ascension  a  and  declination  6  for  the  true  equator  and  equinox  of  date 
can  be  computed  from 


a  =  qatn(y««,x««),  and 

«  =  un-‘  =  fcm-‘  ,  («) 

The  Greenwich  mean  sidereal  time,  Greenwich  apparent  sidereal  time,  Greenwich 
hour  angle  and  local  hour  angle  are  computed  from  formulas  given  on  pages  B3  and  B4 
of  Ref.  2.  The  altitude  and  aaimuth  are  computed  from  Equs.  (10)  and  (11),  respectively. 
The  equation  for  refraction  is  Equ.  (3),  page  B15  of  Ref.  2.  The  equations  for  planetary 
magnitude  are  given  on  pg.  315  of  Ref.  4  with  the  correction  for  Saturn’s  rings  given  on 
pages  362  to  365.  Formulas  for  the  table  on  page  365  were  obtained  using  least  squares. 
Equations  for  the  semidiameter  of  the  sun  and  moon  are  found  on  page  B16  of  Ref.  2.  The 
lunar  parallax  in  altitude  formula  is  on  page  B16  of  Ref.  2  and  the  lunar  phase  formula  is 
on  page  311  of  Ref.  4.  The  lunar  age  approximation  developed  by  the  author  is  within  ±1 
day. 


B.  Sample  ProUam.  Calculate  the  position  of  the  sun,  moon  and  navigational  plan¬ 
ets  for  24  May  1984  at  17^S8™45*  snln  time  at  latitude  38*35^24^  north  and  longitude 
78*18'0(y'  west.  Assume  that  the  temperature  and  pressure  are  the  default  values  of  lO^C. 
and  1010  mb.,  respectively.  The  output  shown  in  this  sample  problem  was  computed  usmg 
Microsoft  BASIOA/D  (double  precision)  mi  an  IBM  PO.  Several  intermediate  results,  such 
as  the  true  position,  and  appaimit  positions  for  the  true  equator  and  for  the  true  equinox 
of  date  are  printed  only  as  a  debugging  aid  for  those  wishing  to  reproduce  the  results  on 
a  computer  other  than  those  for  which  NAVEPBM  is  available. 

Input  Seresoit 

TZAR  (4  DIGITS)  7  1084 

MONTH  NUMBER  7  5 

OAT  OF  THE  MONTH  7  24 
ZULU  TIME  (HH.MMSS)  7  17.6846 
LAT  (DD.MMSS)  (•>-N/-S)7  38.3524 
LON  (DDD.MNSS)  76.18 

TEMP  (DEG  (ZLSIUS)  7  10 
PRESSURE  (MILLIBARS)  7  1010 

PRESS  C  TO  CONTINUE 

The  input  screen  is  used  for  input  only — there  are  no  options  from  which  to  choose.  The 
temperature  and  pressure  are  used  only  for  the  atmospheric  refraction  computation. 
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lat  Ontpat  ScrMmt 


TRUE  POS,  THUS  tCL  k  EQHX  OF  DATE  (UT) : 

63.73236366163271  0 

APP  POS.  TSUE  Ea  k  EQNX  OF  DATE  (UT) : 

63.72674638386701  0 

APP  POS.  TEUE  EQU  k  EQMX  OF  DATE  (UT) : 

4.11446000613666  20.80027008414662 

OREEmriCa  hour  angle  a  DEa  (UT)  : 

00^20.0'  20«64.0' 

ALTITUDE  -  68.60634784607200  68*30.3' 

Aznorra  -  218.6611413863376  218*30.7' 

SD  -  16.8' 

REFRACnOM  -  .4' 

PRESS  C  TO  (2)NTINUE 

The  true  and  apparent  positions  for  the  tme  ecliptic  and  equinox  of  date  printed  are  the 
ecliptic  longitade  and  latitude,  respectively,  in  decimal  degrees.  The  apparent  position 
for  the  true  equator  and  equinox  of  date  is  the  rig^t  ascension  in  decimal  hours  and  the 
declination  in  decimal  minutes.  The  Greenwich  hour  angle  and  declination  are  given  in 
degrees,  minutes,  and  tenths  of  minute  notation.  The  altitude  and  asimuth  are  given  in 
both  decimal  degrees  and  degrees,  minutes,  and  tenths  of  minute  notation.  SD  is  the  sun’s 
semi-diameter.  The  refraction  is  for  the  computed  altitude  of  the  sun. 


Sad  Oatpnt  Scr«enx 
MOOH 

TTttJE  P08.  TRUE  ECL  k  EQNX  OF  DATE  (UT) : 

356.6033313649804  -4.001806604201667 

APP  P08.  TRUE  Ea  k  EQNX  OF  DATE  (UT) : 

366 . 6033313540804  -4 . 901806604201667 

APP  P08.  TRUE  EQU  k  EQNX  OF  DATE  (UT) : 

23.86024633620027  -6.201011880374428 

GREENWICH  HOUR  ANGLE  k  DSa  (UT) : 

164®09.7*  -  6n7.6' 

ALTITUDE  «  5.46172646204884  5^27.1' 

AZIMUTR  -  257.4666216802323  25r28.0' 

PHASE  .31  AGE  -  24  DATS 
SD  -  14. 8^  SD  AUG  -  14. S' 

P  IN  A  -  64' 

REFRACnON  -  10.2' 

PRESS  C  TO  CONTINUE 

The  output  for  the  moon  is  similar  to  that  for  the  sun.  In  addition,  the  moon's  phase, 
approximate  age,  augmented  semi-diameter  (topocentric  semi-  diameter)  and  parallax  in 
altitude  P  IN  A  are  output.  The  true  position  and  the  apparent  position  for  the  true  ecliptic 
and  equinox  of  date  for  the  moon  are  identical  because  no  aberration  correction  is  made 
for  the  moon  (see  Section  A). 


Srd  Output  Sereant 
VKMUS 


TRUE  POS,  TRUE  Sa  ft  EQHE  OF  DATE  (UT)  : 

57 . 72463662418&41  ~ . 6528408164408445 

APP  POS.  TRUE  Ea  ft  EQHX  07  DATE  (UT) : 

57 . 71230545404640  - . 6532006744515710 

APP  POS.  TRUE  EQU  ft  EQNX  OF  DATE  (UT) : 
3.706665485886183  10.01586622158487 

GREENWKS  HOUR  AHGLE  ft  DEa  (UT) : 

06«36.(y  10«01.0* 

ALTITUDE  -  63.67703088674400  63^40.6' 

AZIRUTR  -  227.7063425077327  22r42.4' 

MAGNITUDE  -  -3.4 
REFRACTION  -  .5' 

PRESS  C  TO  CONTINUE 


The  ontpnt  for  Veniu  is  similar  to  that  for  the  snn  except  that  the  apparent  magnitnde  is 
output  instead  of  the  8emi<diameter. 


4ih  Oatpat  Sereonx 
MABS 


TRUE  POS,  THUS  Ea  *  EQNX  OF  DATE  (UT) : 

226 . 1649067006666  - . 7643627638263081 

APP  POS.  TEUE  SCL  A  SQMX  OF  DATE  OTT) : 

226 . 1661434601826  - . 7642106001416870 

APP  POS.  TEUE  EQU  k  EQNX  OF  DATE  (UT) : 

14.80744688277316  -17.30602620686820 

GREENVKa  HOUR  ANGLE  k  DECL  (UT) : 

288*44.3^  -  17*23. 8^ 

ALTITUDE  -  -64.68434644133616  -  64*41.1' 

AZIMUTE  -  62.30703743162624  62*18.4' 

MAGNITUDE  -  -1.6 
REFRACTION  -  -.7' 

PRESS  C  TO  CONTINUE 

The  output  for  Man  is  similar  to  that  for  Vams.  Note  that  the  altitude  of  Man  is  negative, 
that  is,  Man  is  below  the  horison.  The  leads  to  the  anomolous  computation  of  a  negative 
value  for  the  refraction.  When  a  planet  is  below  the  horison,  the  value  computed  for 
refraction  is  meaningless. 
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Stk  Ontpiit  Scracnt 

JUPI71R 

TRUE  P08.  TRUE  ECL  k  EQHX  OF  DATE  CUT) : 

282.0081844M9666  .1413184643412081 

APP  PQ8.  TRUE  KL  k  EQUE  OF  DATE  OTT) : 

382 . 0100408670675  . 1412662061186586 

APP  PQ8.  TRUE  SqU  k  EQNX  OF  DATE  (UT) : 

18.86041325466823  -22.76882700181082 

GREEHWIGH  HOUR  AH6LE  k  J3ECL  (UT)  : 

220«00.6'  -  22«46.6' 

ALTITUDE  -  -61.07026716007840  -  61^58. 2^ 

AZIEUIR  -  206.4711031451716  206*28.3' 

MAGHITUDE  -  -3.1 
REFRACTION  -  -.6' 

PRESS  C  TO  CONTINUE 

The  ontpnt  for  Jupiter  i«  similar  to  that  for  Venus.  See  the  comments  for  Mars  regarding 
a  negative  refraction  value. 


Oth  Output  Screont 

SATTnur 

TSUK  PCS.  TRUE  KL  k  EQNX  OF  DATE  (UD : 

221.6236884013804  2.S80801361264001 

APP  PCS.  TRUE  ECL  k  EQNX  OF  DATE  (UT) : 

221 . 6271566166040  2 . 68001611205868 

APP  POS.  TRUE  EQU  k  EQNX  OF  DATE  (UT) : 

14 . 66067626686222  - 12 . 83630063737880 

GREENWICH  HOUR  ANGLE  k  DICL  (UT) : 

202*17. 4'  -  12*60.2^ 

ALTITUDE  -  -40.04270037701162  -  40*02.6^ 

AZIXUTR  -  60.03733261413306  60*56.2^ 

KAGHirUDE  -  .4 
REFRACnON  -  -.8^ 

PRESS  C  TO  CONTINUE 

The  output  for  Saturn  is  aimilar  to  that  for  Ventu.  See  the  comments  for  Mars  regarding 
a  negative  refraction  valne. 


C.  Program  Lkting— NAVEPHM 


10  BOI  YAMFLASDEUI  ft  PULXKIBBI  PLANET  EPBOCEftIN.  07-20-86. 

IE?  08-27-85  •  0800 

12  8EM  "ASTBO  EPBEN  ?• 

13  lEN  "IIAYEFBN.BA8* 

20  DIF0BLA-N,Q-Y:D1M  A(28),B$(6) 

30  P2-2«An(l) : P1-P2<»P2 :  TP-PI<»PI : ID-PI/180: 8>IO/3eOO 
40  BE-.01720300808:ia$-*1010*:DO«10«:DilMBl(26) 

60  OBPma}-a-ZITCX))*TP:OBFnO(X)-X-TP«ll)TCZ/TP): 

DEFFna)-IITa*M0^ .  6)/M0 
60  DEFnDia)-Z-NO«IITa/NO) 

70  B|(0}-*8DX* :B8(l)-*M0aN* :B|(2)-«?END8« :B8(3)-«KAI8« 

80  B8(4)-«JDPIT11*:B8(S)-«8AT0I>> 

00  GOTO  2000 
100  ; 

110  lEM  2-Aafi  ATAN  PCOi.  IA0IAN8. 

120  A>AI1f(T/a-lE-00*a-0)})-PI*a<0)«TP*CZ>0)*(T<0) tIETUIH 
130  : 

140  UN  AIC81N  ft  AlOCOa  FCIN’S 
180  A8-ATNa/(8qft(l-Z«Z)-lZ-0e«CABa(Z)-l))}:IBI0IM 
160  AC-ATN(8qi(l-Z*Z}/a-lZ-00*(X-0)))-PI«a<0) ;IZnU 
170  : 

180  lEM  HEUQ  SPHEE-TO-IECT  ♦  MOTION  CONPOIATION 
100  lPl/100000:CS-ft*C08(B) 

200  ZB-CS*C08a}  :T1«C8*8IH(L} :  ZB-I*8IN(B) 

210  IF  N-ITIEN  lETOII 

220  10I^1:ZF  NoOTBEN  MO-l«l/IX(K) 

230  PA-B*8qBCAN(N)*lQn/l 

240  Zp-FA*8IN01(H))  ;I^A*8qB(l-K(N)*2)*G08(M(N)):Z-0 
280  A— AP(N):G08UB  280:A— IN(N)  :GQ8DB  200:A^AN(N)  iQOSOB  280 
260  aH-X:VI-T:«B-Z:IETDlN 
270  : 

280  CA<08(A}:8A-8IN(A}:T-X*a-»Y*8A:T^T*GA-X*8A:ZNiT:IETDBN  ;BIN  Z-AZI8  ROT 
200  CA-G08(A}:aA-8IN(A):>T*a-»Z«8A:Z-Z«a-T*8A:T>T:BEI0BN  : UN  X-ACS  ROT 
300  : 

310  UN  RECT-T0-8Pni 
320  G080B  120 

330  L-A2/R0:R>X*X«Y*Y;B-ATN(Z/8qR(R2))/RD:R^qi(R2-»Z«Z)  :RET0IN 
340  : 

380  REN  DECIMAL  TO  ODD  NM.F 
360  ?$-•  MF  X<0TH!»  ?$-•-■  :X—X 

370  X-X-H/1200;T-INT(D:V>-?8+RICHT8(»  ■+8TRia).3)-^*®* 

380  X-600*(Z-Y):Y^IMTa):X$-8TU(1000’fY) 

380  Y8-?8->^MID8(Z$,3.2}-»-."-»IZGBT|(X$.l):REIUIB 
400  : 

410  REN  D0D.NM88  TO  DECIMAL 

420  n-0:F0l  >1T0  LEN(V|}:C|-NIS8(V|.Z.1}:IF  Cl-*. "THEN  IZ-Z 

430  NEXT  :IF  H-OTHBI  X-VAL(?I) : RETURN 

440  X-VAL(LEFT$(V$,n));8N-l:IF  X<0THBI  8N— 8N:X— X 

460  V$-V8-»*0000*:T-?AL(MID8(?8.n<i>1.2})  :Z-?AL(NID|(Y$.  IXt3.2)) 

460  X-8N>»((Z/60-»Y)/60«X):  RETURN 


480  REN  FOHDRNEHTAL  ARCTMENTg 

400  A(i)- .  806434f . 03680110130*78: A(a}-.374807* . 03620164700*78 
500  A(3)-. 260001*. 0367481062*78:A(6}-0:A(7)-. 770073* . 00273700031*78 
610  A(8}-. 003126*. 0027377788*78 

620  A(12)-. 606408*4. 480468670000003B-03*78:A(13)-.  140023*. 00446036173*78 

630  A(14)- . 202408* . 00448040017*78 : A(18)* . 087383* . 00148676328*78 

640  A(16)- .  063866* . 00146661327*78 :A(17)-. 840604* . 00146660466*78 

660  A(18}-.080608*.00023080803*78:A(10)«.086631*.00033060803*78 

660  AC20)-. 814704*. 00023080803«78:A(21)-. 133206*. 00000204371*78 

670  A(2a}- . 882087* . 00000204371*78 : A(23}- . 821218* . 00000204371*78 

680  A(28)-. 400680*3.2604380000000011-08*78 

600  A(4)-A<1)-A(7):A(6)-A<1)-A(3) 

600  m  Z-170  35:A(Z)-FlfAa(Z));lim  Z:RETniN 
610  : 

620  PRIN7  61(0):  REM  80M(0) 

630  X-A(8} :T-A(13) : Z>A(10) 

640  >(6010-17*7J)*8IM(X}*72*8Ill(2*X)-7*C08(Z-Z)*6*8Zlia(l)-A(7)) 

660  L«L*6*8Ilf  (4*X-8*A(16)*3*Z)  -6*C08(3*a-7)  )  -4*81Na-y) 

660  >L*4*Ca8(4*X-8*A(  16)  *3*Z)*3*8IN(2*a-T)  }  -3*8111  (Z)  -3*818(2*  (X-Z)  ) 

670  L-L*8I>*A(7):B-0 

680  R^100014-1676*Ga8(X)-14*C08(2*X) 

600  AM(0)-1 . 00000023: K(0)- . 016700114 :N(0)-A(8) 

700  AP(0}-( (4 . 62778E-04*7J*1 . 710178) *7J*101 . 2208333) *RD: 

IH(0)-0:AN(0)-180*RD 
710  8M(0)-l/320300 

720  00808  100:RBPt:Z8-ZH:78-TH:Z8-ZH:U8-aH::V8-VH:V8-«H:RE7DRH 
730  PRI87  B|(1):REN  MOQK(l) 

740  «-A(2) :X-A(3) :T«A(4) : Z-A(8) 

760  L-33640*8IN(W)-4686*8IN(N-2*y)*2370*8ZH(2*Y)*760*8Z8C2*W) 

760  >L-668*8IH(Z)-412*8IH(2*X)-212*8I8(2*CW-T))-206*8ZH(W-2*Y*Z) 

770  L-L*102*8I8(V*2*T)*166*8I8(2*Y-Z)*148*8INCV-Z)- 128*818(7) 

780  L-L-il0*8IS(W*Z)-66*8Z8(2*(X-Y))-46*8I8(W*2*X)*40*818(Y-2*X) 

700  L-L-38*8I8(Y-4*7)*36*8I8(3*V)-31*8I8(2*W-4*y)*28*8I8(V-2*y-Z) 

800  L-L-24*8I8(2*Y*Z)*10*8I8(8-Y)*18*8I8(Y*Z)*16*8I8(W*2*Y-Z) 

810  I.-L*14*8I8(2* (Y*y) ) * 14*8Z8(4*Y) - 13*818 (3*W-2*Y) 

820  AG-«*16*A(7) - 18*A(12) : L-L- 1 1*818(AG)*0*C08(AG) *4*7J* (COS (AC) *818 (AG) ) 
830  >L*10*8I8(2*«-Z)*0*818(Y-2*X-2*Y)-0*8I8(2*W-2*Y*Z) 

840  L»L-8*8I8(W*Y) *8*818(2*7-2*3) -8*8I8(2*W*Z)-7*8IB(2*Z) 

860  L-L-7*8I8(W-2*Y*2*Z)*7*818(A(8))-6*8I8(W-2*X*2*Y) 

860  L-L-6*8I8(3*X*2*Y)-4*8I8(¥-4*Y*Z)-4*8I8(2*W*2*X) 

870  L-L*3*(8I8(W-3*Y)-8I8(8*2*Y*Z)-8I8(2*W-4*Y*Z)*8I8(W-2*Z) 
♦8I8(W-2*Y-2*Z)) 

880  L-L*2* (SI8(W*4*Y) -8I8(2*W-2*Y-Z) -8I8(2*X-2*y*Z) ) 

800  L-L*2*(8I8(4*W)*8I8(4*Y-Z)*8I8(2*W-y)) 

000  B-18461*8Z8(X)*1010*8I8(V*X)*1000*8Z8(Y-X)-624*8I8(X-2*Y) 

010  B-B-100*8I8(W-X-2*Y)-187*8IR(W*X-2*Y)*117*8IR(X+2*Y) 

020  B-B*62*8I8(2*Y*X)*33*8I8(V-X*2*Y)*32*8Z8(2*W-X)-30*8I8(X-2*Y*Z) 

030  B»B-18*8I8(2*ir*X-2*Y)*lS*8I8(W*X*2*Y)*12*SI8(X-2*Y-Z) 

040  B>B-0*8I8(Y-X-2*Y*Z)-8*8I8(X*A(6))*8*8I8(Xt2*Y-Z) 

060  B-B*7*(-8I8(8*X-2*Y*Z)*8I8(W*X-Z)-8I8(8*X-4*Y)) 

060  >B*6* ( -818 (X*Z) -8X8 (3*X) *8I8(V-X-Z) ) 

070  B-B*B*(-818(X*Y)-8I8(W*X*Z)-8I8(¥-X*Z)*8I8(X-Z)*8I8(X-Y)) 


g«0  M»4*(8Dr(3*V«Z}-8Iira-4*Y))^3*(-8Iiar-X-4*r)-»flIH0r-3«Z}} 
goo  M+3«(-8IH(2*W-Z-4*T)-81]I<3*X-a*Y)«8IS(2*«-Z«3*T)) 

1000  ^B*a*(BIlf(V‘I*2*r-Z)*8m2*W~I-2*r)*8n(3*W~X)) 

1010  BFOOa6208-327746*Ca800-67g04*G08(W-2*T)-46857*C08(2*Y) 

1020  lPE-8g04*C08(2*W)t3866*CQ8(2*V-2«Y)-3337*C08(2*T-Z} 

1030  BFE-2688*C08(Y«2*T)-33S8*C08(W-a*Y<»Z)-2030*CQ8(¥-Z) 

1040  EPE‘»1710*C08(T)-*-ie71*C08(Y-»Z)<»1247«C08(V-2*X)4>704*C08(Z) 

1080  lpE«€2g*C08(2*T<»Z)-6a4*GQ8(V-4*T)4-3g8*C08CV-2*Y-Z}-3e6«C08(3*W) 

1060  EFE-2g6*C08(2*W-4*Y)-263*C08(Y-»Z)f24g*C08(3*«-2*Y)  -2ai*C08(V^2«Y-Z) 

1070  BFE+188*C08(2*X-2*Y)-161*C08(2*Y-2*Z)-»147*C0aaN2*X-2*Y)-142*C08(4*Y) 
1080  B^E’i>130«C08(2*«-2*Y+Z)-118*G08(V-4*Y<i>Z)-116*G08(2«W«2*Y) 
logo  BfE-110*€08(2*«-Z) 

1100  L-L*81HA(1):B-B*80:BfE/23464.8:G08DB  igO:81HL:SETUHM 
1110  ramT  B8(2):BB(  78808(2) 

1120  Z-A(8)  :>AC13) :  >A(14) 

1130  L»(2814-20*TJ)*8IHa)-181*8IM(2*Z)+12*8ni(2*Y)-10*C08(2*a-Y)) 

1140  L-L-»7«Ca8(3*(Z-Y)) 

1 ISO  B-1221S*8IH(Z) *83* (8IN(Y4Z) 48IX(Y-Z) ) 

1160  8^7233S-4g3*CQ8(Y) 

1170  L«L*8D«A(12}:B-B*80 

1180  AM(2)- . 7233316000000001 : K(2)-6 . 773041000000001E-03 : M(2)-A( 13) 

1100  AP(2)-(  (> .  001386380*TJ+ .  60818611) *TJ-»^54 . 3841861)  *ED:  I]((2)-3 . 3g46«ED 
1200  AM(2)-(  ( .  00041*TJ-» .  80088)  *T>78 . 77064722000001)  *  ED :  RN(2)-408S23 .  S 
1210  G08I1B  ig0:BP-E:8ETDB8 
1220  PEI8T  B8(3):EZa(  MAB8(3) 

1230  Z"A(8) : Y^A(16) : Z-A(17) : Q-AC 10) 

1240  L-(384Sl«a7*TJ)*8IX(Y)-»(2238t4*TJ)*8IN(2*Y)‘M81«8IN(3*Y)-62*8Z8(2*Z) 
1280  L-L-22*C08(Y-2*q)-10*8I8a-0)+l7*C08(Y-q)+17*8IS(4*Y)-16*C08(2*(Y-0)) 
1260  L-Ltl3*C08tt-2*Y)-10*8I8(Y-2*Z)-10*8IH(Y+2*Z)+7*Cn8(X-Y)-7*C08(2*X-3*Y) 
1270  L*L-8*8IHa(13)-3*Y)-S*aiH(X-Y)-8*8I8(X-2*Y)-4*C08(2*X-4*Y)*4*C08(q) 
1280  L«Lt3*CQ8a(  13) -3*7)  43*818(2*  (Y-q)) 

1200  Bii6603*8IM(Z)4622*8I8(Y<Z)4616*8I8a4Z)464*8I8(2*Y4Z) 

1300  lPl83031-14170*GQ8(Y)-660*G08(2*Y)-47*Ca8(3*Y) 

1310  L-L«8D4A(18):B-B*80 

1320  AM(3)-1.62368830:BC(3)>.00340488700000002;8(3)-A(16) 

1330  AP(3)-(  (1 . 3128E-04*TJ41 . 06076667) *T>288 . 4317610000001)  *BD : 
I8(3)-1.8407*ED 

1340  AH(3)-( (- 1 . 38gE-06*TJ4 . 7700016700000001) *TJ448 . 78644167) *S0 : 
B)l(3)-3008710 

1380  GQSOB  100:aP-E:BETDBN 
1360  PEI8T  B8(4):8E2fJI]PITEE(4) 

1370  X'-A(18):Y-A(ig):Z-A(a2) 

1380  >ig034*8IlI(Y)46023*TJ4281l4l003*G08(2*Y>6*Z)t601*8I8(2*Y) 

1300  L*L-47g*8I8(2*Y-6*Z)-186*8I8(2*Y-a*Z)4l37*8I8(3*Y-6*Z)-131*8I8(Y-2*Z) 
1400  >Lt70*C08(Y-Z)-78*C08(2*Y-2*Z)-74*TJ*C08(Y)468*TJ*8IN(Y) 

*66*008 (2*Y-3*Z) 

1410  L-L463*C08(3*Y-5*Z)463*CQ8(Y>S*Z)440*8IN(2*Y>3*Z)-43*TJ*8I8(2*Y-5*Z) 
1420  L»L-37*C08(Y)428*8IM(a*X)*2S*8IH(3*Y)'a3*8IJI(Y-6*Z)-18*TJ*C08(2*Y-6*Z) 
1430  L»L*17*C08(2*Y-4*Z)4l7*CQ8(3*Y-3*Z)-14*8IN(Y-Z)-13*8I8(3*Y-4*Z) 
-O*C08(2*X) 

1440  >L40*C08(Z)-0*8IN(Z)-0*8IN(3*Y-2*Z)40*8IN(4*Y-6*Z)40 
•8IM(2*Y-6*Z43*A(2S)) 

1460  >L-8*COa(4*Y-10*Z)47*C08(3*Y-4*Z)-7*C08(Y-3*Z)-7*8I8(4*Y-10*Z) 
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1460  L«L-7*8IN(T-3*Z}'»6*C08(4*Y-6*Z)-6*8Zll(3*T-3*Z}i>6*CQ8(2*Z) 
-4*81X(4*Y-4*Z) 

1470  L»L-4*C08(3*Z}«4*C08(2*Y-Z)-4*G08(3*Y-a*Z}*4*TJ*CQ8(2*Y)<»3«TJ«8ZN(2«Y) 
1480  L»L-»3*CQ8(6*Z)«3*Ca8(8*Y-10*Z)-»a«8ZH(2«Z)-2*8DI(2*X-Y)-»2«8III(2*Z'»Y) 
1400  L-L-2*TJ*81N(3*Y-5*Z}-2*TJ*81N(Y-5«Z} 

1600  ]l^4e02*C08(Y)<»25e*8IB(Y)+227-227*C08(2*Y)^30*TJ*8INCY)-»21*TJ*C08(Y) 
1510  B-fr»16*8IN(3*Y-6*Z3-13*8IH(Y-S*Z)-12*CQ8(3*Y)4>12*8ZN(2*Y) 

1620  ^B«7*G08(3*Y-5*Z)-6*C08a-6«Z} 

1630  BF620883-26122*G08(7)-604*C08(2*Y}«260*C08(2*(Y-Z))-170*G08(3*y-6*Z) 
1640  M-106*8ZHC2*(Y-Z))-01«TJ*8INCn-84*TJ*C08(Y)-»60*8III(2«Y-3*Z) 

1660  lFE-67*8IM(Y-6*Z)'»66*8IN(3*Y-6*Z)-»63*8IN(Y-Z)-6i*C08(2*Y-3«Z) 

1660  W-46*8IMa)-20*C08(Y-6*Z)+27*C08a-2*Z)-22«C08(3*Y)-21*8IlI(2*Y-6*Z) 
1570  I.*L*8D+1:M*80 

1680  AM(4}-6 . 202561 : BC(4)- . 0484042510 : N(4)-A( 10) 

1600  AP(4}-(  (7 . 0406E-04*TJ-» .  6004316700000001) «TJ-»^273 . 2775417)  *60 : 

1N(4)-1. 3031*10 

1600  AH(4)-( (3 . 52222E-04*TJ-<-l . 01063) *TJ*00. 44338611) *RD : BN(4)-1047 . 355 
1610  608DB  100:BP-E:BETDlIf 
1620  PEIHT  B|(6):aEI(  8ATUB1I(6) 

1630  V-ACIO) :Z-A(21) :Y«A(22) : Z>A(26) 

1640  L-2a046*8IN(Y)*6014«TJ-2680*C08(2*V-5*r)*2607*1177*8ZK(2*V-S*r) 

1660  >L-826*C08(2*W-4*Y)*802*8IH(2*Y)-*426*8ZN0r-2*Y)  >220*TJ*C08(Y) 

1660  >L-163*CQ8(2*W-6*Y)*142*TJ*8IN(7)-114*CQ8(Y)*101*TJ*8ZH(2*«-6*y) 

1670  L-L-70*G08(2*Z)-»67*8IiI(2*Z)-»66*8Ili(2*W-6*Y)-^60*TJ*C08(2*«-6*Y) 

1680  L-L*41*8IH(W-3*Y)*30*8IHC3*Y)+31*8I»(W-y)*31*8IlI(2*(W-y)) 

1600  >L<20*C08(2*«-3*Y)  *28*8III(2*«-6*Y*3*Z)i‘28*C08(V-3*Y) 

1700  L-L*22*TJ*8DI(2*W-4*Y)-22*8Ill(Y-3*Z)+20*8IH(2*i-3*Y) 

1710  L-L*20*C08(4*W-10*Y)+10*CQ8<2*Y-3*Z)*16*811I(4*W-10*Y) 

1720  L»L-17*TJ*C08(2*Y)-16*C08(Y-3*Z)-12*8IM(2*W-4*Y)+12*C08(W) 

1730  L»L-12*8Ill(2*(Y-Z))-ll*TJ*8IJI(2*Y)-ll*C08(a*W-7*Y) 

1740  L-L*10*8IN(2*Y-3*Z)+10*C08(2*<W-Y))-»0*8ni(4*W-0*Y) 

1750  L-L-8*8IH(Y-2*Z)-8*C08(2*X+Y)+8*C08(2*X-Y)*8*C08(Y-Z) 

1780  L-L-8*8IN(2*X-T)*7*8IH(2*X+Y) -7*C08(W-2*Y) -7*C08(2*y) 

1770  L-L-6*TJ*8IM(4*W-10*T)+8*TJ*C08(4*W-10*Y)+8*TJ*8I»(2*W-8*T) 

1780  L-L-6*8IN(3*W-7*Y) -6*C08(3* (W-Y) ) -6*C08(2* (Y-Z) ) *6*8I1I(3*W-4*Y) 

1700  L»L+6*8IM(2*W-7*Y)+4*8IH(3*(W-Y))+4*8I1I(3*W-6*Y) 

1800  L-L*4*TJ*C08(W-2*Y)+3*TJ*C08(2*W-4*Y)->3*C08(2*W-6*Y*3*Z) 

1810  >L-3*TJ*8IN(2*X)’»3*TJ*C08(2*W-6*Y)-3*TJ*C08(2*X) 

1820  L-L*3*C08 (3*W-7*Y) +3*0)8 (4*W-0*Y) ♦3*8IN(3*W-8*Y) 

1830  L»L+3*8IH(2*W-Y)+3*8IB(W-4*Y)+2*C08(3* (Y-Z) ) 

1840  L»L+2*TJ*8IN(W-2*Y)+2*8IN(4*Y)-2*C08(3*W-4*Y) 

I860  L»L-2*a]8(2*W-Y)-a*8IH(2*W-7*Y+3*Z)+2*C08(¥-4*Y)+2*C08(4*W-ll*Y) 

1860  L-L-2*8D)(Y-Z) 

1870  B-8207«8IN(Y)  -3346*0)8(Y)  +462*818 (2*Y)  - 180*0)8  (2*Y) + 186 
1880  B-B+70*TJ*C08(Y)-71*O)8(2*W-4*Y)+46*8IN(2*«-6*Y) 

1800  B-B-46*CQ8(2*W-6*Y)+20*8I8(3*Y)-20*O)8(2*W-3*Y) 

1000  B-B+18*TJ*8Z8(Y)-14*O]8(2*W-5*Y)-ll*O}8(3*Y)-10*TJ 
1010  B»B+0*8I8(W-3*Y)+8*8IH(W-Y)-6*8I8(2*W-3*Y)+6*8IH(2*W-7*Y) 

1020  >B-6*0)8(2*W-7*Y)+4*8I8(a*W-5*Y)-4*TJ*8Z8(2*Y) 

1030  B-B-3*0)8 (W-Y) +3*0)8 (W-3*Y) +3*  TJ*8IN(2*W-4*Y) 

1040  B-B+3*ai8(¥-2*Y)+2*8IM(4*Y)-2*Ca8(2*(W-Y)) 

1060  B^65774-63262*0)8(Y)-1878*8IN(2*W-4*Y)-1482*0)8(2*Y) 

1060  Bpa+817*8IN(W-Y)-B30*O)8(W-2*Y)-624*TJ*8IN(Y) 


1070  M«340«8ZH(2*W-S*T)'»347*8ZN(2*W-««Y)<»328*TJ*C08(’n 
1980  lFt-228*8ZN(T}-»140*C08(2*V-8*r}  -126*C08(2*  (Y-T)  } 

1900  M«104«C08(Y-T)-»101*GQ8(2*V-8*T)<»98*CQ8(V-3*Y) 

2000  »>E-73*G08(2*W-3*Y) -82*GQ8(3*Y)«42*8IH(2*Y-3*Z) 

2010  »>t-»41«8IX(2*  (Y-Y)  )  -40«811I(Y-3*Y)  •»40*C08(2*Y-4*Y) 

2020  M-28*YJ-23*81H(Y)i-20*8IH(2*Y-7*Y) 

2030  >L*81HZ:M*8D:1.^L:B0-B 

2040  iM(S)-9 . SS4747000000001 : BC(S)- . 06SS4609270000001 : N(B} -A(22) 

2060  AP(5)-(C9.78642E-04*T>1.08S22060)*TJ-»338.3077722)*RD:ZH(5)-2. 4886*60 
2080  AN(6)-( C>1 . 62181E-04*YJ* . 8731961400000001) *TJ*112 . 7903889} *60 : 
BM(6)-3408.8 

2070  G08DB  100:8P-6:BETD6Y 
2080 

2090  CL8  ;  INPOT  ■YEA6(4  DIGITS)  M 
2110  INFOT  ■MONTH  NDNBEH  "iK 
2130  INPOT  "DAY  OF  THE  MONTH  ■:! 

2160  INPOT  ‘ZOUI  TIME  (HH.MISS)  ■:V|:G080B  420:0T-X 
2170  INPOT  'UT  (00.18(88)  (*N/-8)  ■:V8:GQ80B  420: L>X 
2190  INPOT  'LON  (D00.NN88)  (>Y/-E)  ■:V8:G080B  420:LG-X 
2210  INPOT  "TEMP  (DEG  CELSIUS)  ■.DC 
2230  INPUT  ■PHE8808E  (MILUBAES)  ■:NB 
2240  Q08QB  2720 

2260  J0-3e7*K-INT(7*(X*INTC(M*9)/12))/4)*INT(27S*M/9}*I*i721013. 6 
2260  CLS 

2270  TXJ0-2418020)/36626:0T^1 .04*TJ-16:T^J0*24S1646*DT/3600/24 
2280  T0-T8/30628:T8-T8*0T/24:TJ«T8/3662S-*1 
2290  1M$-*(0T):«:IF  DTHITHPI  T1I8-*(TDT):* 

2300  G080B  490 

2310  G)H<2. 68622E-08*T0*2400. 061336) *70*6.697374660000001*1 .002737909*07 
2320  NO-24 :(a(-FNM(GM):G8-GN-.00029*8INa(5)) 

2330  DL-(*17 . 23- . 03*TJ) *8IN(A(6}}-1 . 27*8IN(2*A(7))* . 21*8IN(2*A(6} } 

2340  DL-0L-.2*8IN(2*A(l}):0L-0L/36OO*6D:BEN  DL-NOTATION  OF  LONGITUDE 
2360  OE-((.00181*TJ-.0060)*TJ-46.846)*TJ*84428.28*9.21*C08(A(5)) 

2360  0E-(0E*.S62*G08(2*A(7}})*8D 

2370  F06  N-OTD  6:QN  N*lGOSaB  620.730,1110,1220,1360,1620 

2380  IF  N-006  N-ITHEN  X-XH:Y-YH:2-ZB:IF  H-ITHEN  0^:y-0;W^):GaTQ  2410 

2390  IF  IMTHEN  O-OH :  T-VH :  IHrH  :X8-X:Y8-Y:Z8-Z:  08-0  :V8-Y:V8-V:  GOTO  2410 

2400  X-X8*XH :  Y^*YH :  >Z8*ZB :  U-U8*UH :  V-V8*VH :  V-V8*VH 

2410  A— IM.:G08UB  280 

2420  P6INT  DN|;*TB1IE  P08.  TIOE  BCL  A  E9NZ  OF  DATE  ■:TM8:G08UB  320:P6INT  L;B 

2430  GD-6:L1-L:B1-B:IF  K-OTHEN  L>L 

2440  (>173.142:BC-6/C:BO(  VELOCm  OF  LIGHT 

2460  X-X-BC*0 : Y-T-BC*V : Z-Z-BC*Y 

2460  P6INT  :P6INT  'APP  P08,  TIDE  ECL  ft  EQNI  OF  DATE  ■;TN$: 

G080B  320:P6INT  L:B 

2470  A— OE.-GOBOB  300:G08DB  320:LP-L:BP-B:IF  N-OTHEN  LS-L:B8-B 
2480  L-L/16 

2490  P6INT  :P6INT  -APP  P08,  TBOE  Eqo  ft  EQNI  OF  DATE  ■:TN8:P6INT  L;B 
2600  M0-360:GH-FNM(1S*(G8-L)) 

2510  P6INT  :P6INT  "GBEENYICH  H006  ANGLE  ft  DBCL  •.TH$ 

2620  X«GH:G080B  360:LS-V$:X-B:G080B  360;P6INT  LS**'':SPC(3):  VI*"'” 

2630  LH-(GH-LG)*RD:CL-COS(LB) :8L-8IN(LH} :BD-B*6D:SB-8IN(BD) :CD-C0S(B0) 
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2S40  L>LT«ID:a-C08(I.D):&MIH(LD):X««8BKI*CB*CL:G080B  160: 
Aft-MiAL-AS/BD 

2560  >-C8*flL:MB*a-0B*6Z«GL:0080B  120:i2N12/BD 

2560  P8IHT  :PEIHT  ■ALniODI  -  ML;  :Z-AL:G08DB  300:PU1IT  8PC(5) 

2666  PBIHT  "AZIXDTH  -  ■:AZ;  :Z-AZ:G08aB  360:PtIHT  8PC(6}:T|W- 

2670  mi  HBOSOB  2780.2780.2780,2860.2870 

2680  IF  PUST  •PHAaS*;!;*  AflS  -"'.Aa:"  DATS 

2500  IF  N>lT]mr  )f>10:P8IHT  ■NAanUDE  *  ■:FinOfS} 

2600  IF  B-0THE6  8-16.004/E:N0-10:PEIHT  *80  -•;8TB$(FHRCS})-^*'" 

2610  IF  HOITHEI  2660 

2620  >23464. 8*t:8-036.75/D;M0-10:PIIET  -80 -■:STB$(FMR(S) 

2630  8-8*(l<»8INCA0)/D):PEIHT  •  80  AQO  -•;8TI|(F)n(8) 

2640  Z>G08ao)/D:G08DB  150:PA-A8/ED*60 
2660  PEUrr  ■?  IB  A  ■  •;STM(FHE(PA))+“'* 

2660  RF-l/TAH((ALi>7.31/CALt44.4})*BD) 

2670  IF^8F*(0O-80)/030)/(li>.0000e«(BF‘^30)*(DC-10)) 

2680  BF-8F- .  00«8IB(  (14 .7*IF-»13)  *80)  :M0-10:  BF^PB8(BF) 

2600  P&INT  ‘BETEACnmi  -  ■:8TB|(BF)-i-*'" 

2700  00800  2720; CL8  .BEH  B.GOTO  2000 
2710  : 

2720  PEIBT  :P8IBT  ‘PBESa  C  TO  CONTINOE* 

2730  d-IBZETl  :1F  C8***TIEB  2730 
2740  IF  Ct-^C^THEB  8ET0UI 
2760  GOTO  2730 
2770  : 

2780  CD-8IN(B8*B0) *8IB(BP*R0) *C08(B8*RO)*C08(BP*ID) *C08( CL8-LP) *80) 

2700  Z«C0:G08QB  160:T*8B*8IM(AC) :Z-00--8K*GD:GQ8DB  120;PA-AB8(Aa/80) 

2600  QB  NGOTQ  2810,2640,2860 

2610  K«(l*Ca8(A3))/2:8B-8IB((L0-Ll)*80):Z-8qi(Z):G0S0B  160 

2820  AG-2B.63*A8/PI:IF  3B>OTEEB  AG-20. 63>Aa 

2630  M>100:1-FBE(Z)  :)I0-1  :AG-FBEa6)  :BETDBH 

2840  NG-(PA*PA*4 . 247E<07*> .  01322)  *PA'*2 . 1715*L06(GD*BP)  -4 :  BETUBN 

2860  MG-.01486*PA-*2.1715*L0Q(6D*BP)-1.3;8ET0RH 

2880  NG-2. 1715*106(60*8?) -8. 03 '.BETUBB 

2870  B0-(168. 1176*1. 304001*TJ)*80:II-(28.0743-.0127001*TJ)*8D 

2880  NN-( ( . 242202*TJ*3 . 08600) *TJ*126. 3620) *80 

2800  J J-( ( . 0171666*7 J- . 4640142) *TJ*6 . 012036) *80 

2000  80-(-.230002*TJ-2. 731270) *TJ*42. 02030 

2010  CB-G08(BO) :Y-8IH(II)*8IM(B0)iC08(II)*CB*8IN(L0-BQ) 

2020  Z-CB*C08(L0-Ba) :G08UB  120:UP«A2/B0 

2030  >BP*80:OC08(0)  :8B-8IN(D)  :8J-8IB(JJ)  :CJ-C08(JJ) 

2040  AG«LP*80*NH:8H-8IH(AG) 

2060  Y«8J«8B*CJ*CB*SN:Z-CB*C0S(AG):G080B  120;l]U-A2/BO 
2070  BB-8J*CB*8B-CJ*8B 

2080  MG-2 . 1716*L0G(GD*BP) -8.68*. 044*AB8(UP*80-DD) -2 . 6*AB8 (BB) *1 . 25*BB*BB 
2000  BETUBN 
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1620-2070 


2090-2240 


Usage 

P2  =  ir/2,  PI  =s  T,  TP  =  2ir,  RD  =  degree-to-radian  conversion. 

EK  s  Sun’s  gravitational  constant.  =  atmospheric  pressure  in  millibar’s 
(1010  is  the  de&ult  value),  DCS  s  temperature  in  degrees  Celsius  (10  is  the 
default  value). 

FNA  converts  angles  in  revolutions  to  degrees  between  0  and  2)r,  FNR  rounds 
to  the  nearest  MO. 

FNM  is  the  X  modulus  MO  function. 

B$  is  an  array  of  planet  names. 

GOTO  20S0  to  commence  execution. 

Two  argument  arctangent  function  with  output  interval  of  0  to  2ir.  The  -lE-9 
must  be  replaced  by  -i-lB-9  in  the  PC-isooA. 

The  arcsine  function  can  be  replaced  inline  by  the  ASN  function,  m  the  PC- 
ISOOA. 

The  arccosine  function  can  be  replaced  inline  by  the  ACS  function  in  the  PC- 
1500A. 

Unscale  the  radius  vector  by  a  factor  of  l£5.  Convert  spherical  coordinates  to 
rectangular  coordinates. 

Return  if  the  body  is  the  moon. 

MU  is  the  reduced  mass  of  the  body.  Compute  the  approximate  heliocentric 
velocity  vector  of  the  body. 

Z-axis  rotation. 

X-axis  rotation. 

Convert  rectangular  coordinates  to  spherical  coordinates.  L  =  longitude,  B  - 
latitude  and  R  =  distance. 

Convert  decimal  degrees  to  degrees,  minutes  and  tenth  minute  notation. 

Convert  DDD.MMSS  or  HH.MMSS  format  to  decimaL 

Compute  the  fundamental  planetary  arguments  [Ref.  10]. 

Compute  the  geocentric  spherical  position  and  vriocity  of  the  Sun.  Convert  to 
geocentric  rectangular. 

Compute  the  geocentric  spherical  and  rectangular  coordinates  of  the  Moon. 

Compute  the  heliocentric  spherical  position  and  velocity  of  Venus.  Convert  to 
heliocentric  rectangular. 

Compute  the  heliocentric  spherical  position  and  velocity  of  Mars.  Convert  to 
heliocentric  rectangular. 

Compute  the  heliocentric  spherical  position  and  velocity  of  Jupiter.  Convert 
to  heliocentric  rectangular. 

Compute  the  heliocentric  spherical  position  and  velocity  of  Saturn.  Convert  to 
heliocentric  rectangular.  Saturn’s  bngitude  and  latitude  are  saved  as  L9  and 
B9  on  line  2030. 

Input  prompts  should  be  designed  so  that  default  values  or  previously  entered 
values  are  displayed.  Pressing  the  return  key  will  input  a  displayed  value.  If 


any  valne  is  changed,  then  the  eniire  valne  most  be  keyed  m.  This  feature 
had  to  be  removed  on  the  TRS-80  ModeM  version  since  it  was  difficult  to 
implement. 

2250  Computation  of  the  Julian  Day  Number. 

2270  TJ  =  number  of  Julian  centuries  from  noon,  January  0,  1900  to  midnight  of 
mput  date.  DT  =  AT  correction  factor.  [NOTE:  This  equation  should  be 
examined  for  accuracy  yearly — see  text].  TS  =  number  of  Julian  days  from 
noon,  January  1,  2000  to  0000  hours  Universal  Time  of  the  input  date. 

2280  TO  =s  number  of  Julian  centuries  from  noon,  January  1,  2000  to  0000  hours 
Universal  Time  (UT)  of  the  input  date.  TS  s  nuniber  of  Julian  days  from 
noon,  January  1,  2000  to  the  input  date  and  time.  TJ  =  number  of  Julian 
centuries  from  noon,  January  0, 1900  to  the  input  date  and  time. 

2290  U  the  AT  correction  is  set  to  aero,  then  the  positions  are  referenced  to  TDT 
(terrestial  dynamic  time)  rather  than  to  UT. 

2300  Compute  fundamental  arguments. 

2310  GM  —  Greenwich  Mean  Sidereal  Time. 

2320  GS  =  Greenwich  Apparent  Sidereal  Time. 

2330-2340  DL  =  nutation  of  longitude. 

2350-2360  OE  =  obliquity  of  the  ecliptic  corrected  for  nutation. 

2370  FOR  loop  to  cycle  throng  Sun,  Moon  and  planets.  N  =  body  number.  The 
loop  en^  on  line  2700. 

2380  X,  Y  and  Z  are  the  geocentric  coordinates  of  the  Sun  or  the  Moon.  Set  the 
velocity  components  of  the  Moon  equal  to  zero. 

2390  Save  the  position  (XS,  YS  k  ZS)  and  velocity  (US,  VS  &  WS)  components  of 
the  Sun. 

2400  Compute  the  geocentric  position  and  velocity  components  of  the  Nth  planet. 

2410  Correct  position  for  nutation. 

2420  Display  true  position  for  the  true  ecliptic  and  equinox  of  date. 

2430  Save  true  distance  GD,  longitude  LI,  latitude  B1  and  solar  longitude  LO. 

2440-2450  Correct  position  for  aberration  (light  time). 

2460  Display  apparent  position  for  the  true  ecliptic  and  equinox  of  date. 

2470  Convert  ediptic  coordinates  to  equatorial  coordinate.  Save  the  right  ascension 
LP  and  declination  BP  of  the  body.  For  the  Sun,  save  right  ascension  and 
decimation  as  LS  and  BS. 

2480  Convert  right  ascension  from  degrees  to  hours. 

2490  Display  apparent  position  for  the  true  equator  and  equinox  of  date. 

2500  GB[  =  Greenwich  Hour  angle. 

2510-2520  Display  Greenwich  Hour  angle  and  decimation. 

2530-2565  Compute  and  display  the  altitude  and  azimuth.  The  altitude  is  saved  as  A9  on 
line  2540. 


2570  Subroutine  call  for  N  =  1  through  5. 

2580  Display  the  Moon’s  phase  and  age. 

2590  Display  the  planet’s  magnitude. 

2600  Compute  and  display  the  Sun’s  semidiameter  [Ref.  2,  pg.  B16|. 


2610  Tyansfer  if  the  body  is  not  the  Moon. 

2620  Compate  and  display  the  Moon’s  semidiameter  [Eel  2,  pg.  B16]. 

2630  Compute  and  display  the  Moon’s  augmented  semidiameter  [Eef.  2,  pg.  B16]. 
2640-2650  Compute  and  display  the  Moon’s  parallax  in  altitude  [Eef.  2,  pg.  B16]. 
2660-2690  Compute  and  display  the  refraction  correction  at  the  body’s  altitude  [Eef.  2, 
pp.  B14-B15]. 

2700  Gimtinue  prompt.  End  of  FOR-NBXT  loop  on  body  numbo-  N. 

2720-2760  Continue  prompt  query. 

2780-2790  Compute  the  phaM  angle  PA.  CD  is  the  cosine  of  the  elongation  [Eef.  4, 
pg.  312). 

2800  Transfer  for  the  MooUf  Venus  or  Mars. 

2810  Compute  the  Moon’s  phase  angle  K  [Eel  4,  pg.  311). 

2820  AG  is  an  approximation  to  the  Moon’s  age  develop^  by  the  author. 

2830  Eound  the  Moon’s  age  to  the  nearest  day  and  return. 

2840  MG  is  the  magnitude  of  Venus  [Eef.  4,  pg.  314). 

2860  MG  is  the  magnitude  of  Mars  [Eef.  4,  pg.  314). 

2860  MG  is  the  magnitude  of  Jupita  [Eef.  4,  pg.  314). 

2870-2980  Compute  the  magnitude  MG  of  Saturn.  This  computation  is  complicated  by 
the  varying  aspect  of  Saturn’s  rings.  Details  foUow. 

2870-2900  Computation  of  the  ring’s  node  BO,  inclination  11,  right  ascension 

NN,  inclination  to  the  mean  equator  JJ  and  arc  SO  from  NN  to  BO.  These 
equations  result  from  curve  fitting  the  table  in  Eeference  4,  page  365. 
2910-2920  Computation  of  UP  s  U'  [Eel  4,  pg.  364). 

2930-2970  Computation  of  UU  =  U  and  BB  ^  ^  B  [Eef.  4,  pg.  345). 

2980  MG  is  the  magnitude  of  Saturn  [Eef.  4,  pg.  314). 
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APPENDECi  The  QATN  Fnnction 


Thia  lontine  is  the  standard  arctangent  fnnction  corrected  for  quadrant.  The  quadrant  arc* 
tangent  fnnction  is  occasionally  implemented  as  the  ATAN2  fnnction,  the  ANGLE  function 
or  the  Rectangnlar-to-Polar  function. 

Entering  variables  are  the  z>  and  y^coordinates,  X  and  Y.  The  exiting  variable 
the  angle  0,  where  — ir  <  0  <  x.  Use  of  the  quadrant  arctangent  Action  is  denoted 
0  =  qatn(y,Z). 

1.  If  X  0,  go  to  step  4. 

2.  Set  0  =  (ir/2)  *8gn(y’). 

3.  Go  to  step  8. 

4.  Set  0  =  arctan{y/X). 

5.  If  JT  >  0,  go  to  step  8. 

6.  Set  0  =  0  +  IT  *  sgn(y). 

7.  If  y  =  0,  set  0  =  x. 

8.  Return. 

Note: 

If  y  >  0  then  8gn(y)  =  +1. 

If  y  =  0  then  a^{Y)  =  0. 

Ify  <0thensgn(y)  =  -l. 

Users  of  Microsoft  BASIC  can  simplfy  the  qatn  fnnction  significantly  by  using  the  code 
^ven  below.  To  return  an  angle  of  0  (designated  by  A  in  the  code)  in  the  range  of  (-t,  t), 
use: 

PI  -  4*ATJI(1):  TP  -  PI  +  PI:  EPS  -  lE-33 
A  >  ATN(Y/(X-EPS*(X-0)))  -  PI*(X<0)*(SGN(y)  -  (Y-O)) 

To  return  a  value  of  A  in  the  range  of  (0, 2x),  use: 

PI  -  4*ATN(1):  TP  -  PI  +  PI:  EPS  -  lE-33 
A  -  ATN(Y/(X-EPS*(X-0)))  -  PI*(X<0)  +  TP*(X  >=  0)*(Y<0) 
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